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_ The behavior of multi-story buildings responding to strong earthquake ground — 


motion has been computed by a high -speed digital computer. _ Response curves 


_and tables are presented, and the effects of certain structural parameters on > 


INTRODUCTION 

multi- building in vibration is an extremely complicated dynamic 
system. It has infinitely many degrees ‘of freedom and is far from linear. 
Many of its components are erratic in their behavior, even in carefully con- 
trolled laboratory specimens. Finding precisely the mot tion of such a system © 
except by actual measurement is beyond accomplishment. . The system can i 
simplified by means of reasonable approximations to reduce it to a solvable 
‘mathematical model. _ However, in doing this, bear in mind that the solution, 

matter how it may be mathematically, is the response of 

“model. he 
Inan the ground shakes irregularly and sometimes violently in 
_all directions. This induces structural response in six components of motion, 
three translational and three rotational. 2 This study is limited toa single © 


a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings: of the 
‘American Society of Civil Engineers, Vol, . 87, No. EM 2, April, 1961, ce. we 
Assoc, Prof. of Civ, Engrg., Univ. of Michigan, Ann Arbor, Mich, 
_ 2 «Methods for Calculating the Earthquake Response of Shear Buildings,” by Robert 
S, Ayre, Proceedings, First World Conf, on Earthquake Engrg., , Berkeley, Calif,, ; June, 
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of structural response, A ‘single translational ‘component of motion may 
induce all six components of structural motion. Even if rocking of the struc- 
tare is prevented or ignored, rotation about a vertical axis will occur unless 
- structure is dynamically symmetric about a plane parallel to the ground 

7 _ motion or about two perpendicular planes. If the structure has a planeof sym- a . 
metry parallel to the direction of ground motion, and if rocking is prevented, | 
_ the response will have a single translational jonmponaat. This is the case con- 


The mass of an actual structure is distributed over ‘the entire structure, 

most cases, , the actual system « can be approximated reasonably well by assum- 
ing that its entire mass is concentrated at several discrete points, Thus a 
—_— mass system is obtained, If the system consists of n point masses, it” 


a the one- »-directional translation considered herein, A multi- story building docs 
_ in fact have mostof its mass situated at the floor levels. _ It can be considered 
a lumped mass system, taking the mass of F each f floor and half the mass of the 
walls and columns inthe stories immediately above and below asa point mass” 
malls at floor level. The error introduced in this approximation can be ex- 
pected to be small in the lower modes of vibration, 3 which are the modes | all 
primary importance in structural res} response to earthquake, 


5 ural resistance of the frame, shear resistance of the > walls, and: viscous opted 
ing resistance. It was assumed that the joints of the rectangular building frame 
were rigid, that is, there could be no relative rotation between the ends of the — 
meer ata joint. In reality, this assumption is not strictly true except per- 
- haps ina frame in which each member at a joint is welded to develop the full 
moment capacity of the member, It was further assumed that all members in 
a the frame, columns, and girders: alike, had ideal elastic- plastic moment- rota- 


Be process for solving the differential equations of motion for such a struc. 


‘The | in a multi- _story. building may have a ‘sizable influence on the 3 re- 

_ sponse of the building to an earthquake, However, the exact behavior of walls | 

is difficult to evaluate, In experimental \ work on determining the periods of 
vibration of the fifteen- story Alexander Building in San Francisco, Calif., John 

A. Blume, F. ASCE, found that the tile partitions offered little or no contribu- 

_ tion to the rigidity of the building at small amplitudes of vibration.9, 8 The 


- brickwork participated in the rigidity, but not to the extent that its modulus of in 


3 Qn the Fundamental Frequencies of Vibration of Rigid Frames, ” by E, F, Masur 2 
_ Proceedings, First Midwestern Conf, on Solid Mechanics, Urbana, Ill., April, 1953, 

4Dynamic Analysis of Elasto-Plastic Structures,” by G, V. Berg and DaDeppo, 
Proceedings, ASCE, Vol. 86, No. EM 2, April, 1960, a 
_ 5 “Period Determinations and Other Earthquake Studies of a Fifteen- ieee Building,” 
John A. Blume, Proceedings, . First World Conf, on Earthquake > Engrg., Berkeley, 
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‘STRUCTURES RESPONSE 
elasticity would imply. A greater participation might be expected at theg great-— 
er amplitudes which would be experienced in a strong motion earthquake, — ry 
r Gw ‘Housner, M, ASC E, has pointed out that a large share of the effect of 
-masonry walls on large amplitude vibrations is because of the energy dissi- : 
_ pated by the masonry units cracking and grinding together. 7 Even plaster par- 
titions may be able to dissipate appreciable amounts of energy. 
EX eriments conducted at Stanford University on brick and concrete shear 
; walls® suggest that the load-deflection characteristics of such walls can be 
7 reasonably described by a tri-linear elastic- semielastic-plastic curve, Tests a 
involving strain reversal in masonry shear walls have not been reported, to : 
: the writer’s knowledge. For reasons of mathematical convenience, an elastic- - 
plastic type load-deflection curve, with hysteresis on reversal of strain, was 
used to represent the overall behavior of a shear wall, This is admittedly. a 
‘rough approximation. No claim is made that masonry has. elastic- — prop- 
ertiesinshear, 
“snr The evaluation of the shear resistance of the walls was seeiitias by an 
incremental technique similar to that used for the elastic- plastic frame, - Be- 
cause of the absence of carry-over from one story to the next, the details of 
applying the technique were for shear walls. 
STRUCTURES ANALYZ ED 
computations \ were made for buildings of two different. heights, 
nd for two different | earthquake components, The buildings chosen for the © 
omputations represent steel framed structures proportioned for the lateral | 
- wind loads specified intypical United States building codes. Framing diagrams © 
are shown in Figs. 1 and 2. The f frames were not designed for: seismic forces, 7 


a 
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the four-story building, were made using four different 
mathematical models to represent the structure. The first of these was an un- 

- damped elastic system in which the resistance wastaken asthat dueto theelas- _ 
tic frame alone, with no allowance for damping forces, inelastic frame defor- 
mation, or the resistance of walls and partitions. Second was a damped elastic 
_ model in which a resistance due to inter-floor viscous damping was added to 
the elastic resistance. In the third ‘model, the system was taken as the un- 
‘damped elastic resistance of the bare frame, plus an elastic- plastic type com- 
ponent of resistance due to shear walls, For the fourth model, the system was _ 

7 takenas anelastic-plastic bare frame, without shear walls or viscous damping. 
For each of these systems, the elastic resistance of the bare frame was the | 

dominant resistance and in the last three cases it was modified by 

additional ‘component. For the eight-story building, the analyses were 


limited to the elastic undamped case and the elastic frame with interfloor vis- 


The mathematical models used for the response computations were based» 
on frames one bay wide, equivalent to the three- bay frames in the sense that | 


the total flexural stiffness and ultimate flexural strength of all the columns in 


7 a 7 “Behavior of Structures Dur ing Earthquakes,” by G. W. Housner, Proceedings, Vol. ;> 
Behavior of Reinforced Concrete Shear Walls,” by J. Benjamin and Ay will- 
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April, 1961 
each story were the same in the equivalent frame as in the original frame. 
The same equivalence applied to the beams in each floor. The reason for this - 
substitution was that in the elastic- -plastic case, computing the moments and 
hinge rotations for the ne acon three-bay frames would have required more 


the 2,000- word memory of the ‘computer, Even though the c original three- bay 

frames could have been used for all cases except the elastic-plastic case, the 

equivalent frames were used for all analyses in avoid un- 

‘the different systems. It turned out that the eesiniies frames were slightly 

stiffer frames, the difference being of the andes of 1%. 


FIG, 1.—FRAMING DIAGRAM, FOUR-STORY BEN 


be noted in the aforementioned sense not 


the elastic- -plastic case. For example, elastic- -plastic r multi- bay frame 
it would be possible (indeed, almost certain) that at some stage of the response, 
"part of the columns in a story would be in the elastic range while other col- 
umns in the same story ‘contained. plastic hinges. In the single bay frame, 
_ der the assumptions used herein, plastic hinges occur simultaneously in sym- 
- metric pairs, and all columns ina story reach plasticity at the same time. ¢ 
a was assumed for the viscous damped models that interfloor viscous damp-- 
ing forces of a magnitude of 3 kip sec/in. were effective in each story. This — 
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STRUCTURES RESPONSE 

For the four-story building with shear walls, the same elastic- plastic type 
“resistance function was used for the shear walls in each story. The shear 
stiffness in each story was taken as 30 kips per in, and the ultimate shear | 


- strength in each story was takenas 20 kips. The shear wall stiffness was pur- 


187. 


kips 
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FIG, 2. —FRAMING DIAGRAM, -STOR 


- posely taken to be less than the stiffness of the elastic frame, and the shea 
wall strength was taken to be quite low. This was done to obtain an idea of ho 
_ the response is affected by adding a small amenend of seal wall resistance t 
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earthquake accelerograms were in the analyses; namely, the E-W 
7 _ component of the earthquake recorded in Helena, Mont., on October 31, 1935, | 
a and the N- -S component of the earthquake recorded at El Centro, Calif., on May 
18, 1940. The were assumed to be true records of the actual 


approximated quite closely by piecewise — functions, which are also con- 
- venient for use in the computer, The time-acceleration coordinates of the 
uy intersection points of successive line segments» were put in punched card form 7 
for m machine input data. . Fig. 4 shows an enlarged portion of an accelerogram 
record and a plot of the same portion of the punched card accelerogram. 
4 There is no zero axis on the accelerogram record and the true axis cannot 
be located with certainty. For these studies the axis was selected parallel to 
_ the horizontal baseline appearing on the record and located such that if ‘zero 
- initial velocity were assumed, _ the computed ground velocity at the end of the 


= punched card accelerogram would also be zero. Subsequent investigations of 
wy 


‘been assumed to be very stiff a 
4 would have dominated the response. 
q 

= = = a 


"STRUCTURES RESPONSE 7 
_ computed ground displacements have led to better methods of locating the axis.9 
_ Differences are slight, however, being of the order of 0.003 g. While the ef- 
of such differences on computed ground motion may be be substantial, their 
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COMPUTED RESULTS 


a Ther response of the four~ ~story undamped elastic system to the Helena ground a 
“motion is is shown in Fig. 5. After the first 2.5 sec, the response tends to be 
7 nearly periodic, with a period of 2 0 sec. The natural period for the funda- 


¥ mental mode of this system is in fact 2. 01: sec, The irregularities in the de-— 
flection curves are evidence of vibration in the higher modes. 
7 _ The response of the four- -story viscous damped system to the Helena ‘ground 
- motion is shown in Fig. 6. As might be expected, there is a marked difference _ 
_ between this and the undamped system, Here the deflections buildup to a peak 
at about 2.9 sec, whichis about the end of the most severe portionof the cad 
‘shock. _ Thereafter, the structure oscillates with a period of about 2 sec, and 
_ with steadily decreasing amplitudes, The irregularities in the deflection cur 
" ves are not as pronounced here, indicating that the higher modes of vibration 
Since the viscous damped system can be ‘solv ed by the method of normal — 
_ modes, it is of interest to evaluate the first mode of response and | compare . 
with the total response obtained by direct numerical integration. The fraction 
of critical damping is 0.148 in this mode. . The first mode response is shown > 
in Fig. 7. It compares favorably with the total response shown in Fig. 6. How- | 
ever, , the comparison cannot be expected to be this goodinall cases. © 
Fig. 8 shows a plot of the response to the Helena quake for the four- “story 
-undamped structure with shear walls added, shear walls act as conserva- 
_tive elements when they are in the elastic range, and become dissipative ele- 
“ments when the elastic range exceeded, One effect « of the ‘shear walls ap- 
‘pears as a shortening of the period of vibration of the structure. After the 
"structure reaches its extreme amplitude, it oscillates with a period of 1.45 
sec, or about three quarters of what the period was when the shear walls were 
not present. The fundamental period of the | system is 1. 45 sec when the shear 
walls are inthe elastic range. The effect of the dissipative forces om to out 
the . amplitude does not approach zero as it did inthe case of the viecous dam- 
_ ped system, but rather diminishes to the point where the oscillations have a 
_ amplitude s sufficiently small to permit the e shear walls 1 to remain in the elastic © 
range and, thereafter, the system behaves as an undamped system. The dis-— 
_ sipative forces do not tend to damp out the higher modes of vibration, = 
ie No response is given for the elastic- plastic frame to the Helena ground mo- © 


tion, because the ground motion was not severe enough to induce plasticity in Pt 
Fig. 9 shows the deflection-time curves for the four-story undamped elastic 
structure responding to the El Centro ground motion, The El Centro earth- 
quake was considerably more intense and of longer duration than | the Helena 
earthquake. Only the first 10 sec of response is shown in Fig. 9, The char- 
acteristics of the response are much the same as appeared in Fig. 5 for the 


“ “Integrated Velocity and Displacement of Strong Earthquake Ground Motion,” by 
Berg and G, W., Housner, Seismological Soc. of America be 
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same system and the Helena ground motion, ‘The amplitude of oscillation in- 


=. creases steadily during the time shown, The © computations | were e actually car- 


er than that shown in Fig. 9. The response shown here was computed on the 
7 assumption that the system would remain elastic regardless of the he amplitude — 


of the deflections. This would be the case if the frame members wer were built of. 
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st ip! 


(6) Plot of Punched» Card 
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: an extremely high elastic limit. ASTM- al steel has an average 
yield strength of approximately 42,000 psi under dynamic load conditions, and 


- 4 _ would be strained weil into the inelastic range by the oscillations shown in Fig. 
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; _ For the s same ne framework with elastic-plastic frame ae of 42, 000 | psi 
yield strength, the response to the El Centro ground motion is plotted as de-— 7 
_ flection versus time in Fig. 10. At first glance, these curves appear almost yn 
~ unrelated to the curves for the elastic system, except for the first 2.5 sec of 
response, However, a close inspection reveals a great dealof similarity. The 
_ periodicity of the oscillation is nearly the same as for the elastic system, How- — ; 


ever, in the elastic-plastic system, the oscillation occurs about a position of — 7 
equilibrium which shifts away from the zero position and keeps changing as in- | 
- elastic deformation progresses, Moreover, t! the shift in equilibrium position is — 
not the same for all floors. Also, the amplitude of oscillation about this shift- _ 
ing equilibrium position remains nearly constant, ‘whereas in the elastic sys- > 
Although the total ‘deflections, ‘measured from the zero are “greater 
for the elastic-plastic system than for the elastic system, the amplitudes of 
the oscillation are only about half the elastic amplitudes. Here, as in the ear- 
lier case where shear walls were ¢ considered, the system behaves asa con-— 
servative system ) during the part of its response when the system is elastic, © 
behaves as a dissipative system plasticity is present, The 
 — of plasticity induces a decrease in stiffness of the structure, which 
amounts for the increase in maximum deflection measured from the zero po- 
sition, Plastic deformation also dissipates energy that would otherwise 
stored in ‘the structure as recoverable strain energy, and thus acts to reduce 
_ the amplitude of oscillation, The locations at which plastic hinges formed dur 7 
ing the response are shown in Fig. Lae 
The maximum values obtained for various meen in the analyses are 
gummarised in Table 1, The important parameters are the maximum deflec- 
tions: at the shear coefficient for or each story S 


Maximum shear in in ith story 
Weight of structure 


_ Maximum shear in ith story 
Maximum shear in first — 
the seismic ‘coefficient | for ea each § story or 


Maximum shear in ith story ( 


(3) 
Weight of structure above ith story ie 
_ ‘These e parameters are Shown in a. 1(a) for the Helena ground motion for — 
iy r the four- -story building. Table 1(b) shows the ‘same parameters for the four- 

story building and the El Centro ground motion, and Table 1(c) shows the E 


rameters for the eight- “story and the Helena ground motion. ‘For 


J 
aig 
...2) 
a | 

7 San Francisco building codes are also tabulated, 
oa _ _In the analyses reported in this paper, the displacement response of a dam- oer 
“ ; ped linear system is approximated reasonably well by the response of the first a 
Pe _ mode alone. For each of the analyses, the maximum displacement in the first oe 


STRUCTURES 


mode differs from the total maximum m displacement by about 10% at most, For _ 


_ However, there is no reason to believe that this would be true in general. For 
the same earthquake, the deviation is greater for the taller structure, | The 
modes beyond the first would be expected to be more important for a taller 
structure, simply because there are more of 

_ In terms of seismic coefficients, which are the parameters of the first sig- 

nificance, the first mode approximation is less attractive, especially in the up- 
_ per stories. For the four-story structure, the deviation of the first mode re- _ 
_ Sponse from the total response interms of seismic coefficients is 11% at most. — 
the eight - -story structure, the deviation ranges from practically zero in the 
bottom three stories to 44% in the top story. 

te Any form of energy dissipation in the system tends to reduce the amplitude _ 
of oscillation below that which would prevail if the dissipative forces were not: 
there, This is true even for the case in which the energy dissipation occurs 
through the deformation of plastic hinges in the structural member. The frame 


becomes stiff the plastic hinges” are e present, but the plastic 


as can seen clearly from Table results here partially substan- 


> Housner’s thesis is that if a structure canbe built in such a way that inelastic 

_ deformation progresses far enough to dissipate a substantial amount of energy, _ 
then the structure can be designed on the basis of lower seismic coefficients _ 
than would have to be considered for a comparable elastic structure, This ae 


| 


because the « energy dissipation w — act to reduce the magnitude of the ‘maxi- 


an Information now available (1961), is not sufficient to enable the ven 


mum shears in the structure. _ Asie 


both the structural and n non- aaecmeesd: elements, and i in many buildings sthe 1 non- 
: structural elements are apt to be the greater contributors, It would seem on 
2 the damping forces in the structure, if they a are to be considered as viscous — 
damping forces, might be reasonably well represented by interfloor 
= or symbolically by interfloor yr dashpots. 

It is interesting to note that the same system of interfloor damping forces 

produced 14.8% of critical damping in the first mode of the four-story struc- — 7 
ture, but only 6.25% of critical damping for the eight- -story structure, — 
the damping forces are produced inlarge part by non-structural elements, they — 
might not be significantly different for similar buildings of different amped 
4 If | this is so, it can be expected that tall buildings ar are l less aia damped 
than short buildings of similar type. 

_ The actual damping forces in a structure are, of course, not truly viscous. 
Fort the most part they are attributed to internal friction within the = compon- 

ents of the structure. L. S. Jacobsen has shown10 that :a single degree sys- 

tem with non- viscous damping, subjected to a sinusoidal input force, can be 
_ treated as a a viscous damped system with without ut any | appreciable loss — accuracy, 


10 “Steady Forced Vibration as Influenced by Ds nnping.” by = S. Jac 
tions, ASME, Paper APM- 15,1930, 
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each cycle ¢ of oscillation isthe same for the viscous damping forces as for the 

~ actual damping forces, But in the earthquake problem an additional difficulty 
appears, for the fraction of critical damping for equal energy dissipation must — 

vary as the amplitude of the motion varies. This is easily demonstrated for 

the single degree c oscillatory system by ‘considering the energy dissipated per 

eycle of oscillation, 11 The same tendency exhibits itselfin the effect of the in- 


elastic deformation of the shear walls. When elastic-plastic Shear walls were 


added to the ‘undamped s: system, the reduction in maximum response was ‘much 
more pronounced in the more intense earthquake. If the energy dissipation “7 
4 the action of the shear walls were to be represented as equivalent viscous dam- 


_ intense ground motion, This agrees with the results of large amplitude vibra- 
tion tests reported by T. , Hisada and K, Nakagawa, 12 which indicate that the 
fraction | of critical damping as determined experimentally tends to increase 
_ A good reason for the desirability of reducing a structure to an equivalent 
venom damped system is that this puts the structure in a condition to be a 
alyzed by response spectrum techniques. Recent efforts to extend spectrum 
_ concepts to non-linear systems show great promise, ‘but are not yet — ad 
general use on multi-story structures.13 
The reader is cautioned against ‘drawing too-broad conclusions the 
limited studies presented herein, Some of the computed results might be sub- 
stantially altered if the properties of the structures were changed. a 
tic-plastic case, the single bay frame used in the computations is not dynam- 


case, the solution can be no truer than the mathematical model used in it. _ 


—. 


7 ically equivalent to the three-bay frame from which it was derived. Andi in se 


The primary function of dynamic analysis in earthquake engineering is to _ 
~ lead to a better understanding of response phenomena, to an improved evalua- 


on the response. 
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“The velocity distribution far upstream from the contraction is 
pment by a sine curve. The solution of the governing differential system _ 
the streamline pattern at the immediately 


INTRODUCTION 

_ The separation of flow and formation of eddies at ‘the corners immediately — 
preceding an abrupt contraction in a conduit, in which a homogeneous viscous 
fluid is flowing, can be attributed to the rotationality of the flow far upstream — 
and the subsequent action of viscosity. To demonstrate the effect of upstream 
_rotationality on eddy formation, Yih2 obtained the exact solutions of the 
- governing differential equation for the cases (a) of the two-dimensional flow 
_ with a line sink symmetrically located with respect to the longitudinal bound-— 

: : ary, and (b) of an axis-symmetric flow with a point sink located on the axis of os 
_ symmetry. In these solutions, viscosity was considered only as far as its ef- 

_ fect in producing rotationality of the flow flow is concerned, and it it was neglected in 


= up the differential system. 


_ Note.—Discussion open until September 1, , 1961, To extend the closing date one month, — 
<a a written request must be filed with the Executive Secretary, ASCE, This paper is part - 
of the copy righted Journal of the Engineering sGeohentos Division, _Proceedings of the 
American Society of Civil Engineers, Vol. 87, No. EM 2, April, 1961, aon ee 

_ 1 Engrg. Research Div., Scott Paper Co, , Chester, Pa, 
Solutions for Inviscid Rotational Flow with Corner Eddies,” Cc. s. Yih; 


Journal of Fluid Mechanics, Vol. 5, Part 1, January, 1959, p. 36. oe = => 
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Flow ina two-dimensionalchannel with an asymmetrical abrupt contraction 
d 


pril 1961 
‘The » differential system, and its solution, for the flow in a two-dimensional | 
conduit withan asymmetric contraction are presented herein, The asymmetric 
ay contraction first is approximated by a line sink, and second by a finite slot open- 
_ ing. In both cases the streamline pattern . obtained from the solutionof the ‘dif 
_ ferential system shows the anpemetion zone (corner eddy) at the region im- 
mediately preceding the contraction. A comparison of the two solutions indi-- 
cates that for large contraction ratios, small slot openings compared to chan- 
3 width, the contraction can be approximated by a sink without changing the 


FLOW INTO AN ASYMMETRIC LINE SINK 


‘The flow under consideration is a steady two- -dimensional flow ina a long 
4 channel of width equal to unity and gaping in a wall with a line sink lo- 
Cat cated ed asymmetrically at at the corner. uch an arrangement approximates that 


0 


02 


the flow under a sluice gate, located at the end a with t 
= opening replaced by a line sink. The trace of the line sink is taken as the 
“origin of a Cartesian system of coordinates, with | the y-axis extending along : 


lower boundary of the conduit 1). et mentioned prviasly, at 6, 


&g | 
— 
core 


CORNER 


equation’ governing th the two- ang 


Eq. 1 hows that ‘the | Z- nt of the vorticity al ong 
q. shows tha -componer ) y along a a 
stant and is therefore afunctionofv only, 
If the flow of a viscous fluid is considered, the velocity distribution far up-— 
stream from the sink should be parabolic or that forturbulent flow. But either 
of these distributions would make Eq. 1 non-linear and, consequently, rather 
difficult to obtain a simple solution. Instead the velocity distribution far up 
stream is approximated by a sine distribution which is nearly parabolic, = - 
approximation makes Eq. 1 linear. Such an assumption can be tolerated since 
the effect of non- uniformity of the upstream velocity on the flow pattern in the 
~ vicinity of the line sink is sought. If Umax is the velocity at the center line “7 
“the chan channel, that is, the dimensionless upstream is 


(ty) dy =— cos (1 y) 


The absolute value of P is immaterial and, hence, the additive term - Vn ion 


be neglected. Substituting Eq. 3 in Eq. 1, fora section far ‘upstream, 7 1 


ut cos(m y) = 7 from Eq. us 
which is valid not only far upstream but also everywhere downstream, =— 


the vortici tyr remains constant along a stream line, and viscous forces are neg- 
lected. Eq. 1then becomes = 

(6; 

and is is valid for the \ whole field of flow, far ne as well as) near the line 


ei v is dimensionless, and L is the length ‘a the channel, the boundary c con-_ 


) 


“Hydrodynamics,” ” by Lamb, Dover New York, 1945, p, 244, 


| 
misthen 
5 
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"substitution of Eq. 7 in Eq. 1 


and the boundary conditions tee are 


conditions for were obtained by consideration of Eq. 7 and 


1 


= 


essary inorder to ‘oitein a Fourier series type of solution, ess 9 satisfies Eq. 
. = and the boundary conditions (i), (ii), and (iii) for ¥y. The coefficients Cy 
and C, are determined from the boundary condition (iv) for vy in in the following 


s (1 y)] y) dy = - ( 


L 


‘Thus, the stream function defining t the flow patt in a ‘channel of le 


“Fourier Series and Boundary Valve Churchill, McGraw-Hill 


14, 
and (n 7 y)dy 
e—_ 
— 
Cc. E (n?-1) — 


CORNER EDDIES 

4. ‘The flow pat pattern for an assumed L of 5 is shown in Fig. 1(a), in which ‘the. 
‘Separation of the flow and the formation of the corner eddy are evident . The 


; - point of ‘separation is. at x x =-5, because at this section the velocity has been 


‘assumed to have sine distribution. 
INTO AN ASYMMETRIC FINITE SLOT 
‘The geometry of the flow considered herein is the same as that in the pre-— 


vious case except that the line sink at the origin has been replaced by a two- 
dimensional finite slot opening of widthB. 
‘The velocity far upstream is again assumed to have a sine distribution, while 


the dimensionless horizontal component of the velocity the slot, is is as- 
and vanishes at y= it has Cat y= 0, if is Lagrange’ 
stream function, the value of nil is determined with the use of the continuity. : 


| 


cos 


which the function values are taken at x the previous 
‘section, for a sine distribution of the velocity f far ‘upstream, the value of ¢" is” 


=—cos(my).... 


COs 


q 
ii) 


Again, in order take. care of the conditions near the slot opening, assume 


é 
as) 
asx--L and 0S y < 
= for x <0 
“a 


“The bo boundary « conditions for in thi this case 


and 


and for x = 0 and 0< sys 
the horizontal component of the SSunieies at the slot and the boundary condition _ 


= -Ccos sx) a dy cos(ry). 
The last boundary is 
“ein 
forx=Oand0<y<B, : 
i Solution of Eq. 20 by the method of separation n of variables : yields — 


v= +) sin(r 


sin(n ty). 


‘which is the same as Ea. 9 and s satisfies 20 and ‘the preceding 


The coefficients Cy and are determined from conditions 
and (v) for in the following — 


sv 
in - 1+ cos(r y y)] s 


4 


me 
&g (i) wt = 0 — 
g & 
= 
a 7 L(n -1) sin(n7y y 


CORNER EDDII EDDIE 


‘Thus, the stream nadine, y’ "defining the flow ‘pattern in a channel with a a 
4 cos(7 B) 


The flow pattern for an assumed L of 5 an 
“which the corner eddy is again evident. a 


The analyses presented clearly indicate that | the formation of corner eddies 
immediately preceding - an abrupt contraction can be attributed to the — 
- ality of the flow. Thetheory of irrotational flow® does not predict the existence | 
corner eddies for flows with similar contractions. 


wl A comparison of the two flow patterns in Fig. 1, indicates that for B= 0,15, 
7 bd B — 0, there is small difference in the size of the corner eddies. _ There- - 


charge openings, of width small compared to the channel width, large” 
> ‘discharge openings it is expected that the corner eddy will be smaller, | 
_ The asymmetric location of the line sink, or slot opening, made it necessary © 
- to specify a finite value of x at which the velocity has the assumed distribution, — 
if a simple Fourier series type of solution is desirable. Similar limitations © 
are necessary inthe ° solution of some problems of mathematical physics. This 
- limitation, however, was not necessary in the cases considered by Yih’ be- 
| cause the sink was symmetrically located with respect to the flow. | aa 
The solution of the basic differential equation, Eq. 1, for actual velocity | 
"distributions, that is, when Eq. 1 is non-linear, can be accomplished with 4 
“use of computing machines, or by the method of relaxation. It should be kept 7 
_ in mind that solutions of Eq. 1 do not represent the flow pattern§ inside the 
eddies, because this flow does not originate far upstream. Consequently, these 
os solutions a are likely to be valid only for the regions outside the corner eddies 


and fc for flow flows of viscous fluids at it high F Reynolds numbers, These limitations — 


- 


“Theoretical L. M, Milne- -Thomson, The Macmillan Co., 


(6 “On Steady Lanioas Flow with Closed Streamlines at Large Reynolds Number, . by 


- Batchelor, Journal of Fluid Mechanics, Vol, A, 1956, p. 117, 
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rary, <= is, in producing rotationality of the flow he flow only, 


‘The wi writer expresses his and indebtedness to C. S, 
Yih, of the University of Michigan, for his * numerous — in carrying 
_ out the analyses presented in this ‘note, whe 
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BRITTLE EXPLAINED BY NEGATIVE RESIDUALS 


E. Grinter, 


The ‘hypothesis presented is that residual lateral | compressive stresses 
‘may suppress ductility across a notch and permit cleavage fracture to occur. © - 
_ A theoretical explanation it is given. _ Tests made at Brown University on = 


med to check this They provide substantiation within the range 


— 


‘Ultimate load design of structural steel 
to make a lack c of understanding of brittle fracture a factor of critical impor-_ 

tance, Any clarification of the subject, , particularly by use of the _ simple con- 

- cepts of mechanics used by structural designers, will have a strong influence | 
toward improved design. Hence, the writer has attempted to explain the only 
static tests of brittle fracture that have been widely circulated and discussed © 
through use of the elementary concepts of stress and strain, both wantis a and 

_- The principle established is that lateral compressive residual stresses, in- 

: _ troduced by precompression across a notch, crack or flaw, under load rever- 

_ salmay suppress necessary ductility across the notchand thus permit a Cl a cleav-_ 


Note, .—Discussion open until September 1, 1961, To extend the closing date one month, 7 
a a written request must be filed with the Executive Secretary, , ASCE, This paper is part — 
_of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 
of Civil Vol. 87, No. EM 
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_ Slight plastic flow in the x direction will as readily eliminate residual stresses 
perpendicular ‘thereto - in the y and z directions, However, because this argu- 
_ment is not essential here it will not be emphasized. = = — | 


es Steel will fracture in two. distinctly d different modes, that is, by shear flow 


and by tensile » cleavage. Shear | slip leads s to the ductile 1 mode of fracture, How- 
ever, temperature is also a factor. The 1.0 shear-tension ratio of pure “in 
_ sion may be compared with the 0.5 shear-tension ratio for uniaxial ill 
_ The former permits a ductile fracture at temperatures far below the brittle | 
fracture temperature range for simple tensile specimens. 
{: As the shear-tension ratio reduces, the temperature at which a cleavage 
‘ fracture may occur increases until at a zero shear-tension ratio ( equal tri- 
a axial tensions) only a cleavage fracture is possible. A sufficiently low shear- 
tension ratio to initiate a brittle fracture ‘May occur at a notch because tension 
applied a across a notch induces lateral tensions along a and perpendicular to to 


notch, Fractures always initiate at a notch or at a notch-like flaw. 


PLASTIC FLOW INCREASES TRIAXIAL TENSIONS AT END OF CRACK > 
> «nila A plate may resist tensile loads perpendicular to each other | that produce a 7 
biaxial stress field. At an internal crack or microscopic flaw a condition sim- 
ilar to that shown in Fig. 1 will then exist. The stresses acting on the micro- 
— scopic particles A and B located at the end of the crack or flaw must be con- — 
sidered, This problem has been analyzed by elementary elastic concepts on 
the assumption that the particles A and Bmust remainof the same thickness, 2 
: The - conclusion reached was that there is an induced tension perpendicular to 
the plate at its mid-depth at particle A equal to u Syield, that is, the product 
of Poisson’s ratio and the yield stress of the material when the Oy stress at 
the exact end of the crack in Fig. 1 approaches t the yield stress, No distinc- 
_ tion will be made herein between the numerical values of yield stress and elas-_ 
_ For oy = 37,000 psi and » = 0,3, the particle A in ‘Fig. 1 may resist tensile — 
stresses of 37, ‘000 psi in two directions and 11,100 psi in the third perpendic- 
ular direction. Such a particle will | have its range of possible ductility re- 
"stricted but but not totally ‘suppressed, Tests have > shown that a specimen such 
as Fig. 1 ‘must be subjected to arelatively low temperature before the ductility | 
_ may be sufficiently reduced to permit the flaw or microscopic crack to extend 
mie Close inspection of brittle fractures in large structures confirms the ron 
that some slight plastic flow had taken place before the brittle fracture oc- 
curred, The stress concentration at the end of the crack in the vertical direc- 
tion in Fig. 1 las n oy when n may reach a factor of ten or more, _ When n Oy 


“How Does Propagation of a Small Crack Take Place in a Thin Plate?, L. 


4 rinter, , Welding Journal Research Supplement, 1945, p. 319: 


pyo 
_ It will not be necessary t nce of a simplified stress- ome: 
a strain diagram for structural steelrepresenting perfect elasticity followed by eal 
e ideal plasticity. Nor will it be useful to review the unquestioned effect of the ey 
: _ initiation of plasticity in eliminating residual stresses of eithersign whenthey im 
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2 BRITTLE FRACTURES 


“reaches Ovje ield the: deformation €y €y at the exact endof the crack will grow rela- 

: tively large, Therefore, because plastic flowtakes place without change in vol- 

ume of the plastic material, particle A must reduce its thickness rather sharply 

Z perpendicular to the plate in order to elongate in the vertical or y direction, | 

_ Particle B does not undergo the plastic strain Eye Nevertheless, because B is 
attached to A it must also reduce in thickness while resisting the action of A. . 


as for it evident that plastic deformation at the 


Microscopic Microscop! 


Particle A B 


LS Bk 


Cylindrical 
gerne! _ 
Dimpies on he 
rface of Plat 


FIG. 2. —DRUCKER- MYL ONAS NOTCHED SPECIMEN 


end of the evack: will force particle A to develop tension and B compression, _ 
‘perpendicular to the plate. Because plastic deformation is much greater 
elastic deformation, it is concluded that: the tensile stress developed | perpen- 
dicularly to the plate at the end of a crack (Fig. 1) should considerably ae 
the elastically induced stress of or 100 psi, as soon as plastic 


EM 
FIG, 1.—TRIAXIALITY OF STRESS AT END OF CRACK 
IN PLATE ACTING UNDER BIAXIALLOADING 
{Surface of Plate | | 
g 


aes plate due to the initiation of plastic flow across the | crack may ‘conceivably 
approach Cvield thus producing nearly equal triaxial tensions. would help 
’ Bays the observed maggocantn of ductility that accompanies a running crack 


MECHANICS OF BRITTLE FRACTURE OF 


YIELD TENSION P/A 


The notched plate, Fig. 2(a),is a convenient specimen for testing the typical 
effects of a stress raiser onthe action of a plate under uniaxial loading. . Ten- 
-sile tests of notched plates have repeatedly produced brittle fractures when nthe 

nominal P/A stress has reached or slightly exceeded the virgin yield ‘stress. 
of unnotched material. The test specimen of D. C, Drucker, F, ASCE and C, . 
_ Mylonas4 shown in Fig. 2(a) usually fractured by such “cleavage at -10° F_ 

when P/A reached yield tension, The fractured specimen developed a small _ 

depressed area or dimple on each surface of the plate surrounding the notch, 7 
This surface indication of a cylindrical kernel of plastic ee leads | to an 
understanding of the internal mechanism of embrittlement. 

_ Long before the tensile P/A loading in Fig. 2(a) approached the initial yield 

stress across the throat of the plate between notches, the stress concentration — 

6 at the notch had reached the yield stress. Plastic deformation across the 
“notch must occur without change in volume of the plastic material. Hence, in | 
order for the plastic cylindrical kernel of Fig. 2(b) to increase in diameter in a 

the y direction it should reduce in diameter in the x direction, This lateral L 
‘deformation induces atension 0, at the level of the notch resisted by a 7 

; sive reaction stresses in the x direction above and below the cylindrical ker- 7 
" nel, _ The stress arrows, in Fig. 2(b) indicate action of the notch kernel upon 

the surrounding material, that is, tensile forces pull inward, However, the — 
_ appearance of dimples on the surface of the plate shows that the internal _ 
area of the plate | must be very resistant to flow toward the notch in the x di- a 
rection and that considerable flow inward along the notch occurs in the z di- __ 
The complexity of tensionand shear variation in and around the cylindrical 
kernel of Fig. 2(a) is such that only the crudest computations of the internal 
tensile stress 0, necessary to produce a surface dimple can be made. How- — 
ever, the assumption that shear flow will take place around the surface of the 7 ; 
i cylindrical kernel due to tensile force on its central circular cross-section, — | 

Which is comparable to the initiation of an indentation under the action of a rivet 

punch, provides the relation o7 ,d= =2 T ‘Here, t is the plate thickness, d 
the diameter of the cylindrical kernel or surface dimple, and 7 is the ‘unit 

shear that produces shear flow. For 7 = 7 5 Oz = 0.5 Oyielg, as in plastic flow : 
under uniaxial tension, d= t is obtained. Therefore, the diameter of the dim- . 
ple would need to be equal to the thickness of the plate for a surface dimple 

to be produced by a yield point tension oz at the mid-depth of the plate. Be- 

cause th the e diameter of th the dimple i in a structural plate is usually less than the 


Po 3 “Triaxial Tension at the Head of a Rapidly Running Crack in a Plate,” ” by i, Mi. 
Frankland, Journal of Applied Mechanics, December, 1959, p, 570, 
4 «Static Brittle Fracture Initiation at Net Stress 40% of Yield,” by C. Mylonas, 1 D.C. 
Drucker, and J. D. Brunton, Welding Journal Research Su lement, 473 - 
>» 
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be considered, These were analyzed by J. M. Frankland? 
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"thickness, it appears that the value of 07 not only equals Cyjejq but that it oc- 
curs 's from the mid- depth of the plate outward to within a “x from the sur- 
- face ; about equal to the radius of the dimple. In contrast, the compressive re- -_ 
action stresses surrounding the kernel will be mainly elastic because the cross- 

_ sectional area of reactive material is much larger than the cross-sectional area 
For a thin sheet (in which d > t) the stress 0, would doubtless be far be- | 
low  Oyield: Accordingly, there would bea sufficient difference between the _ 
principal stresses to produce internal shears that should initiate plastic flow 
inthe y direction, Only as the thickness increases to structural dimensions, ; 
above 1/4 in., will the lateral resistance to deformation in the z direction be- 


4 441 ait 


come sufficient for to approach yield resulting possibility of brittle 


fracture. re, 


q 
| | 


NT 
TO PLASTIC PRECOMPRESSION 
‘Through extensive tests Drucker and Mylonas have shown that the 
hotched plate and the simple cylindrical without a notch, when precom- 
pressed plastically a as indicated in Fig. 2(c), n may later fracture by cleavage 
if loaded lightly in tension.5,6 The precompressive deformation in Fig. 3(a) 
was from 1% to 3% measured across the notch, The precompressive deforma- 
tion of the cylinder in Fig. 3(b) was 40% to 50% longitudinally. This was thought — 
to represent the maximum possible unit strain at the root of the notchin Fig. | 
due to stress concentration, Precompressed specimens of both types 
tested in tension at -10° F fractured by cleavage at a low tensile P/A stress, 
However, only a few of the cylinders cracked below the yield tension while a | 
large proportion of the notched and precompressed plates fractured at a P/A 
stress well below yield, Hence, the embrittlement due to precompression was 
muchmore severe in the notched plate specimens than inthe cylindrical speci- 
mens, The residual stress pattern in each type of specimen will be developed 
to determine whether the brittle fractures observed can be explained by the 
- action of negative residual stresses in suppressing plastic deformation, Note 
in Fig. 3 that the direction of arrows represent action of the notch kernel on 


= surrounding material, that is, compressive residual forces act outward and 
tensile forcesinward, 


_ RESIDUAL STRESSES IN NOTCHED PRECOMPRESSED PLATE 


notch by a compressive veforce in they direction preninonert small raised pimples 

on each surface at the notch, _ Accordingly, compressive plastic flow at the 
“notch” took place along both y and z axes, It may be assumed without much 

_ chance of error that some e plastic flow also took place outward from the notch en 
along the x axis. This action is represented in Fig. 3(a) by flattening or in- _ 
creasing the horizontal axis of the imaginary plastic cylindrical kernel at the } 


— Hence, under the compressive loading, triaxial compressive stresses | 
_ 5 On the Exhaustion of Ductility of E- Steel in Tension Following Compressive Pre- 
strain,” by D. C. Drucker, C, Mylonas, and G, Lianis, Report, _ Div. of Engrg., Brown 


- 6 “prestrain, Sizes and Residual Stresses in Static Brittle Fracture Initiation,” by Cc. 
_ Mylonas, Welding Journal Research Suppl lement, 1959, Pe 414-s, 
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~s exist at the mid- point of th the note notch with each axial stress ress equal tothe yields stress” 


4 in compression, — When the compressive load is removed, the cylindrical plastic 
_ kernel will tend to increase its vertical diameter, its horizontal diameter, and 


its length by elastic recovery, 


| - Expansion of the kernel cannot occur freely because it is surrounded by a 
7 much larger volume of elastic material. For example, 2 as the cylindrical ker- 
nel increases in length from outward plastic flow in the z direction its motion 
is resisted by surface shears that in turn develop tension in the surrounding 
material equal in total magnitude to the total compressive stress of the 
inthe z direction, Of course, the average unit tensile stress 0z' developed | 
_ may be only a fraction of Oz, te compressive stress, (Fig. 3(a)) because the 


me 


a7 RESIDUAL STRESSES IN PLASTIC 


“FIG, 3. \L STRESSES LEFT BY BY 


"resistive material surrounding the kernel may have several times its volume, — 


4 Note particularly that the compressive and tensile lateral forces are in dynam- a 
ic equilibrium during plastic flow and will remain in static equilibrium when 

. flow is terminated. Hence, there is little opportunity for release of residual i 
lateral compression in the cylindrical kernel in the z direction, . Although the 
situation inthe x direction is much less readily visualized, there will bea resi- 
— dual compressive stress 0x at the level of the notch that will be in equilibrium 
with residual tensions 0x" above and below the notch, The important conclu-— 

_ sion | is that the residual lateral compressive stresses of yield value required 7 
for lateral plastic flow are essentially frozen into the plate at the level of the 


notch by a condition of stress equilibrium, It i is not how appreciable 


| 
= 


relaxation of stress or occur, Hence, 
the biaxial residual compressions must represent large fractions of the yield : 
stress, probably 0.75 yield or greater. 
oa In the vertical or y direction the entire » plate will recover « elastically « due 
_ to removal of the compressive loading. The plate will elongate in the vertical 
direction. The only resistance to this motion is from the plastically shortened _ 
material across the root of the mre that must be stretched by a residual ten- 
‘sion of considerable magnitude. As required by statics (ZFy = 0) the high 
residual tension | Oy at the notch (Fig. 3(a)) must be balanced by alow residual 
compression oy' across the major part of the throat of the plate between the — 
notches, Hence, the residual stress pattern at the mid-point of the notch, as 
- shown in Fig. 3(a),. is that of compressive stresses equal toat least 0.75 yield : 
_ for Ox and 07 with perhaps an equal tension 0, across the notch, _ 


STRESS ARROWS | 
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FIG, 4 WITH IMAGINARY SLICES ACROSS THE NOTCH ILLUSTRATING 
STRESS ADDITION IN OF LOADING AND LATERAL STRESS © 


bar > 
A study of the precompressed cylinder Fig. 3(b) « establishes the ¢ icine 


hoop tension. on, The re residual stress Oy is is zero, 0. 
_ BIAXIAL RESIDUAL COMPRESSIONS SUPPRESS DUCTILE DEFORMATION 


— sion and the triaxial tensions at the notch produced by a final tensile load ap- 4 
plied across the notch, For * clarity only the stress arrows showing t the action 


of the cylindrical kernel on its etait are ‘included. It can be 
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observed that the residual stress and the stress concentration due to vertical 
load tend to add numerically in the ¥ direction because they are of the same 
sense. Laterally, in either the x ors - direction, the load stress is of opposite | * 
er sense to that of the residual stress, Because the stress concentration factor 
y= for the notch in the y direction in Fig. 4 is several times P/A, and because a 
large residual stress of the same sense and direction already exists due to - 
elastic recovery of the plate after precompressive plastic flow, it is clear that — 
a very low average P/A stress of perhaps 10% of Cyield will be likely to pro-— 7 


duce a resultant yield stress in the y direction at the notch, , At the same time 


sat 


_ the lateral residual opepenamionn and their accompanying shears are reduced, | 


volume of the we material. The cancellationo or r reduction of lateral stresses. a 
and their accompanying shears by opposite residuals suppresses ductility. — wi 
A Perhaps" a better way todescribe this phenomenon is is to note that the lateral 


stresses must change elastically from about -0.75 Syield te to +1. 00 yield 
variation of 1.75 0 yield During the same instant of time the cross- noteh 


_ sion can only increase elastically. If its range of stress is comparable to that — 
of the lateral stresses, it should instantaneously reach 2.5 Oyjejq if the resi- 

dual tension 0, was equal to 0.75 oyjejq. Evidently, the cleavage stress under 
these conditions may be less than the instantaneous elastic tensile stress with 
the result that a brittle fracture is initiated, Whether the fracture crack — 
minates quickly or progresses is on the magnitude of the overall 


In tension tests on notched plates, after an initial p iailaiiindidaitiia: ‘it has 


De oe reported that nearly all specimens develop at least a short brittle frac- 
Ber crack at low load that may be heard even if not observed on the surface, 


Mylonas a short crack at the may then the remainder of 
Plate to stretch The 1 reason, the concept previously 


a _ Fora perfect precompressed cylinder r (Fig. 3(b)) a P/A tensile stress ss equal 
to °yield would have to be applied | to develop a brittle fracture, because there 

_ is no notch stress concentration ni nor any residual tension in the vertical di- 

_ rection, _ However, there are flaws in all materials. An occasional specimen 
will iene a flaw within the central core that would produce a stress concen- _ 
tration both under precompression and later under an applied tensile load. Ac- 

cordingly. in agreement with tests, some cylindrical specimens might be ex- 


pected to fracture at a final P/A well below OW Oyield> it appears: there- 


> 
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7 stresses so that the material at the end of the crack does not have its ductility _ ae 
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BRITTLE FRACTURES 


a= stresses explains adequately the brittle mioegaas observed in a teats of of steel 1 


| 


under precompression followed by load reversal. 


a Notch Sevbvittiement and Structural Fractures. —Tests have shown that 7 
‘notches or other severe stress raisers will produce brittle fractures under an 


: application of asingle tensile stress equal tothe yield P/As stress, For mem- 
bers not subject to stress reversal, having determinable stresses well below 
the elastic limit, notches" should not endanger the structure because their ‘stress 
concentrations will be reduced sharply by the initiation of plastic flow. an 
For material subject to stress reversal, such as ship plates and certain 
members of of bridge trusses, a sharp notch may produce | a high compressive 
— “stress concentration that will upset or thicken the notch material by plastic 
_ deformation even under a normal P/A compressive stress, On release of the 7 


there ‘may remain lateral residual stresses at right angles 


expected to approach the yield stress. Then, ‘whe the load is reversed, the 
lateral tensile ‘stresses must first eliminate the lateral residual compressions 
before. they can induce the shear stresses necessary for plastic flow. The de- 
- day or suppression of plastic flow under tension across the notch even for an © 
- instant may permit thetensile stress at the root of the notch to reach the ulti- 
mate cleavage strength of the steel so that a brittle fracture will result. ~— 
A comparable upsetting action followed by stress reversal can occur from 7 
welding. | The weld heat, by expanding a restrained plate locally, may upset 7 
material at a notch or flaw that may remain below an annealing temperature, - 
Then on cooling, the main force at the notch will ill automatically be reversed to © 
tension with the possibility of brittle fracture, 
a It must be recognized that the phenomena of plastic flow andcleavage frac- an 


ture are subject to many variables within the material itself. It is therefore - es 


expected that geometrically similar test specimens or comparable structures — t 
sien owe considerable variation in method of fracture and inthe fracture load- 
evidence that ductility is being suppressed, 
An inspection 1 of visible notches should make it possible to locate upset or 
thickened material in a structure, In some cases it may be possible to drill 
out such material leaving a circular hole that produces a much lower stress 


at the notch, Doubtless expert use of a torch could be made to release the resi- 
| compressions and thus eliminate the immediate danger of brittle fracture. an 
However, the. existence of the notch remains as a threat if stress reversal can 


a _ Tests are needed to determine whether stress reversals repeated many a 
at a low P/A compression along with a high stress concentration factor 
_ may eventually suppress ductility in the same manner as occurred for the tests 
considered 
@: 
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BENDING OF BEAMS ON ISOTROPIC ELASTIC SOLID 4 


B B. vesié, 1M. ASCE - 


SYNOPSIS 
ae In the first part of this paper, Biot’s solution for bending, under a -concen- ‘: 
trated load, of an infinite beam resting on an isotropic elastic solid is brought © 
to a complete form and extended to the case of loading by a couple. Integrals 


in both solutions are evaluated ane an- 


: of subgrade reaction k. For the ‘first time the realm of validity of the well- 
‘ known Winkler’ s hypothesis, and the amount of error induced by its application, ; 


— It a shown that this hypothesis is s practically satisfied for beams of infinite 
length. Concerning beams of finite length a criterion is given for — 
of cases where the conventional analysis may give reasonably accurate results. 
An expression is derived for evaluation of the coefficient k inappropriate cases. 


= 


Notation. 


: they first appear, in the illustrations or in the text, and are arranged alpha- 


for convenience of reference, in Appendix II. 


_ - Note.—Discussion open until September 1, 1961. To extend the closing date one month 4 7 
a written request must be filed with the Executive Secretary, ASCE, This paperis part | 7 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 


_ American Society of Civil Engineers, Vol, 87, No. EM 2, April, 1961, fame 
Prof, of Civ. Engrg., Georgia Inst, of Tech., Atlanta, 


- 
a 
a 7 ‘The second part of the paper contains a thorough investigation of the relia- | ce 
af the conventional annrnach to the nrohlem means of the ene 
| | 
ia 

| 
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The problem of bending of nan elastic solid ch 
acterized by a Young’s modulus Eg, and a Poisson’s ratio ‘Ls has been treated 
_ by several authors. M. A. Biot2. derived a rigorous solution for a beam of in- in- 
finite length acted upon by a concentrated load. E. E. De Beer, 3,4 A. Habel, 5 
‘ _J. Ohde, 6 De Beer and D. Krsmanovitch,? M. Kany,® and others presented ap- 
bres proximate solutions for particular cases of beams « _ length loaded by « one 


- Biot considered the problem of bending, under a cunconiveted load P ( Fig. 1) 
a flexible beam of a semi-infinite homogeneous, elastic, 


a FIG, ae INF FINITE BEAM ACTED U UPON BY A — LOAD 


isotropic solid. For the bending M at any point x of the 


in which Qa denotes a dimensionless parameter; b ‘is the half- width of the beam; 


to a function of bafe (c = 1.00 if if the distribution of pressures a across the width | au 


2 “Bending of an Infinite Beam on an Elastic Foundation, ” by M. A. Biot, . Journal of — 
Mechanics, Transactions, Amer, Soc. of Mech, _Engrs., ‘Vol, 59, 1937, pp. 


- 3 Le calcul de poutres reposant sur le sol,” by E,E, De Beer, Annales des Travaux 


Publics de Belgique, Vol, 51, 1948, pp, 393-418, 525-552, 653-691, 721-749, 


th Computation of Beams Resting on Soil,” by E, E. De Beer, Proceedings, Second 
- Internatl, Conf, on Soil Mech, and Foundation 1 Engrg., Vol. “1, » Rotterdam, 1948, pp, 119- 


“ NSherungsberechnung des a uf dem elastisch- -isotropen Halbraum aufliegenden 


va elastischen Balkens ,” by A, Habel, Der Bauingenieur, Vol. ‘19, 1938, pp. 76-78, 
_ 6 “Die Berechnung der Sohldruckverteilung unter Grandungskérpern,’’ by J. Ohde, 
Der Bauingenieur, Vol. 23, 1942, pp. 99-107, 122-127, 


7 


“Calcul de poutres reposant sur le sol, ” by E, E, De Beer and D. Krsmanovitch, | me 
naa des Travaux Publics de Belgique, Vol. 52, 1951, pp, 981-1018; Vol. 53, 1952, no 


“Berechnung von ” by M. Kany, Ernst, Berlin, 


Tae 
| 
ia 
| 
— 


| 
SOLID. 


= the beam is uniform; 1. 00 <C <1. 13 if the ‘deflection across the width of the 
beam is uniform); Epl denotes the stiffness of the beam (Ep being Young’s — 
= and I the moment of inertia of the beam); Es, pg are Young’s modulus — 
and Poisson’s ratio of the subgrade respectively; 8 = ba k refers toa dimen- 
sionless parameter; (8) is a function tabulated for 8 > 0.1 and given by an 
assympotic expression for 8 <0.1; and P denotes the concentrated load. Hav-— 
ing computed the approximate value of bending moment under the load, that is, _ 
' 7 for x= 0, he found that, in the interval 0.01 <b/e <i, withan error of the order _ 


| 


= avd 


“will be the same as obtained by the more rigorous theory (Eq. 2). 
“modulus of subgrade reaction” is sometimes used for the coefficient of sub- 
2 grade reaction k which represents the constant ratio of contact pressures | and a 
sesso along the beam, and has a dimension of pounds per cubic inch. ih 7 
= paper the term “ "modulus® is strictly reserved for K=kB, which isa quan-— 


should give finite values in spite of the infinite value of the integrand fora = 0 
Ss However, bending moments other than the maximum one as well as other 
influences in the beam (deflections, shearing forces, contact pressures) were 7 
not evaluated by Biot. Also, no solution forthe case of loading by a couple was — 


‘The results which follow inthe next three sections of this paper are oe 


Young’ s modulus E. ) Biot: also showed that the expression for deflections 


BEAM ACTED UPON BY A CONCENTRATED LOAD 
| Eq. 1 is taken as a basis for complete solution of the considered problem. 7 
Using the well-known relationships between the bending moment M(x ) on the — 
_ one hand and the shearing force > V(x), contact pressure p(x), as well as the : 


— 
4 
_Mmax = 0.332 Pb | Cll - | 
means of the conven- 
—_ 
-=2kb=kB 
a 
— 
% 
| 
= 
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(x 


sin(a x/c) x). 


bing f a cos(a x/c) a costa x/c) da ‘Jol 


a 
2 x ‘c) da. 


a 
+ 
x) © cos (a x/c) 2) da 
Jq(x) =- f x/e) ae = = - {cos cos(a x/c) da + f 
and that the divergent part of the integral can be identified with zero.9 
_ a Eqs. 5 through 9 enable the computation of any of the influences in an in- 


4 finite beam acted upon by any concentrated loads. 


‘UPON BY A COUPLE 


le [ Fig. can be tow by u: using a an approach similar to that used in 


the conventional theory of subgrade reaction.10 | 

Consider a beam ‘[ Fig. 2(b)] under the action of two equal and opposite 


- loads P at a distance Ax. Any influence Z will mae 
> [zy (x + Ox) Zy (x)] 
9% *Grillage Beams on Elastic Foundation,” by B. . Drapkin, Proceedings ASCE, 


on las No. 771, Vol. 81, August, 1955, 
_ 10 “Beams on Elastic Foundations, ” by M. Hetenyi, “Univ, of Michigan Press, Ann 


> 


> 
— 
— 
| 
| 
: 
— 
thy 


SOLID 


in which is the for P= “1, The? loading a couple M may be con- 


sidered» asa case the loading shc in in Fig. ‘provided that 


Thus, multiplying me dividing: the right sic side of Eq. 11 by Ax al and passing to the 


dx 


_M 


w(B) cos (a x/c)d 


a3 + 


These expressions, together with Eqs. 5 9 a complete 
-_ eral solution of the problem of bending of infinite flexible beams ona somi- 


ov 
© EVALUATION ON OF INTEGRALS 
In view application of derived expressions one of the six integrals Jn (x) 
mt be evaluated numerically for different abscissas xf = §. The remaining 
he integrals may then be found by | integration o or differentiation of Jn, provided | ; 
_ their value in one point (such as at x= 0) is also so computed. For ‘complete _ 
evaluation the integral Jo(x) (Eqs. 7 or 14) was chosen. 
_ The numerical integration of this type of integr ral evidently presents er 
+ difficulties | which can be overcome in different ways. In this case an integra- 
tion by interpolation was found to be suitable up to a certain point M ( Fig. 7 


from which point on the integrand c could have | been replaced v with s sufficient a ac- 


i 
— 
| 
. Applying this relation the following expressions are derived from Eqs. 5 to 9: 7 _ 
| 
is) 
p(x) = — = ..... (17) 


1961 
b 


—AN INFINITE SEAM UPON A COUPLE 


HE INTEGRAL Ja) = = 


— 
. 
EVALUATION OF | 


ELASTIC SOLID 


- first part t (0 <a using s 


rule whereas for the: part the following approximation was applied: 
a cos (a x/e)da cos to 


3 | 


The position of the point M as well as . the number of points for interpolation — 
after Simpson were chosen So as to ¢ assure an seeatmtah of i 2% in the final value’ 


TABLE 1,—VALUES OF INTEGRAL (ax/e) 4, da 


+0, "45101 +0 66837 
0.03704 05006 =| + 0,20586 

-0,08960 | -0,05947 
 -0,14779 
- 0.13826 
- 0.05811 
—-=0,01366 


of the integral. The results of the > computation are given in Table ‘4 ant also - 
“shown ny on 4. 


"ANALYTICAL EX EX PRESSIONS FOR INTEGRALS 
," ohady of the shape of curves shown in Fig. 4 revealed their | cian ke 
larity with the curves of ' damped wave type. After some investigations it was - 
- found that the computed apeeumieeal values of the integral Jof fit the analytical 


= ae 
o(x) J9(0) en [cos x x - sin X x] 


4 

_in which ated analytically, so 
0,75  =0,08041 
4  =0,00183 = 
i 


42 1961 
in ‘in which J9(0) represents ; the e value of the integral for x=0 Bie 


a parameter having the dimension (length) and defined 


in which x0 e abscissa of the first zero of ‘the integral. 


FIG, 4.—VALUES OF THE INTEGRAL Jolx) 


_ It was found that ‘the value of this abscissa can be approximated, witha 1 max- 


imum error of the order of (2% in the interval 0.01 <b/e <1, = the expression = 


isa function of b/c | only, it can easily be proven that the. five 
Jo(x) to J5(x) may be found from the expressions: 


= ) (cos x sin x) 


X sin 


EM 
a 
= 
: 


J3(x) =33(0) cos 
- = 3 * (cos x + sin’ x). 


‘mere 310) and the values of the integrals for x = =0 and 
‘spectively. One of these integrals” is already ki known: 
are 


‘The remaining ‘four of them were computed numerical 
by the integrand by its asymptotic values. The com- 


ion o and d Js principle, miler to that of 


2,—NUMERICAL VALUES OF INTEGRALS Jo(0), Jy 4 


4828 


tw was found that the following z analytical expressions approximate the c com-_ 


‘puted values of integrals, so that the maximum error in the interval 0.01 <b 


- <1,00 remains of the order of 4% for Jo(0), 5% for Jy ; 3% for J4(0) and 


0. 333 


.. (2) 
a 
— 
pe 
on : be 
a 
7 


= 0.250( 


Eqs. 21, 22, : and 25 through 35 analytically _—— with sufficient accurary, the 
integrals” appearing in the general solution by Eqs. 9, and 


through 17 of the considered problem. 


RE 


LIABILITY OF THE CONVENTIONAL 
nani and mat foundations, as wellas grillage beams, resting on ne 


:, i behavior of | which is 3 usually well simulated by that of elastic Solids, and 


are based on the assumption that the contact pressure at any point of the ‘beam 


is proportional to the deflection of the beam at the same point (Winkler’s hy- 
_ However, theoretical investigations mentioned earlier in this paper, as well 
as experiments ‘made by De Beer!! and others showed clearly that the pres-— 
_ sure distribution at the contact between beams or slabs and a may be 
quite different from that assumed by the conventional analysis. Winkler’s hy- : 
pothesis | is, , thus, generally 1 not satisfied. And yet Biot has shown that itwas _ 
possible, in solving a particular problem of bending of a beam of infinite length | 
acted upon by a concentrated load, to choose a certain value of the coefficient ap | 


of subgrade reaction k, which will give a bending moment diagram in reason- | AG 


_ able agreement with that obtained by treating the subgrade as an elastic solid. — 


_ Therefore, the second part of this investi are was undertaken with the he aim 


Is Winkler’ ever satisfied beams resting on an elastic- | 


+ jould be nesigned to the coefficient of subgrade reaction k ? ees = 
Pa... . Whi * amount of error is to be expected if the analysis is made by the el- 
cementary theory based on Winkler’s hypothesis? 
4, in. 2 the analysis by the conventional theory furnishes results which in 
‘cases agree reasonably while in n other cases” deviate considerably from 
the results obtained by ¢ a more rigorous analysis, will it not be possible to es- - 


VALIDITY OF WINKLER’S HYPOTHESIS 


To answer the first question it. is ‘sufficient to compare Eqs. 4s. 26 and 27 with | 
Eqs . 29 and 30, respectively. It can be seen that the ratio of the contact pres-_ 7 
sure and the deflection of the beam does not depend on x but on b/c only. _ 
«Tests for the Determination of the Distribution of Soil Reactions Underneath 


Beams Resting on Soil,” by E, E, De Beer, Proceedings » Second Internatl, Conf, on > 
‘Soil Mech, and Foundation Engrg., 2, Rotterdam, 1948, pp. 142-148, 


— 
44 = April, 19610 
a 
must be some cases wereh the coj q 
a 
ae 
i= 
— 


a infinite length vastili on a semi-infinite elastic subgrade. Any problem of 
bending of ; of an infinite beam having a stiffness Ep I and a width B = 2 b and rest- 
ing on a semi-infinite subgrade defined | by a Young’s modulus Eg and a Pois- 
son’s ratio pg can be treated with reasonable accuracy by the conventional an- 


4 
alysis using a a coefficient | of reaction k 


VALUE OF THE COEFFICIENT OF SUBGRADE REACTION Bye. | 


What value should be assigned to the coefficient The ‘conventional 


_ ysis leads to solutions which contain, as it is well known, the same damped- 

wave as Eqs. 21 and 26 through 30, with a factor 


36) 


p-curves obtained by the analysis: will become congruent to the 
_ corresponding curves obtained by the rigorous solution of the same problem. 
biases ‘Eq. 25 to 36 and neglecting the small power of C (1 - Hz) ial 


the maximum bending moments (Eq. 4). Introducing B insteadof 2b and taking 
 C=1.10asa oiticientiy correct value for any practical purpose, the following 


expression for Ka is obtained: 


ea! B= K.= 
ghee 
- From Eqs. 37 and 38 ke can be computed if Eg is known. In some instances © 
the value of the coefficient k for a square plate having a width B: is known, De- 


fae kB=K and knowing that, for a rigid square plate: 


‘the folowing can ‘be 


— 

a : ‘This ex ression differs somewhat from that found by Biot by simply equating i 
= 


tt is interesting to compare this supression with the formulas presented by 


Terzaghi,!2 Hon. M. ASCE and F. Vogt. ** 13 The Terzaghi formula gives 
= = 0.67 K irrespective of the flexibility of the beam. Vogt formula gives Kx 
0; for beams of finite length L it gives” irrespective of the of 

0.421 


the possible range of values of | han 


where the lower limit is is for very rigid be beams and the upper limit is for very 


CONVENTIONAL ANALYSIS 

- ‘ What amount of error is to be expected if beams of infinite length on on an e- 
_ lastic- isotropic subgrade 3 are re analyzed | by means of the elementary theory based 
_ Using the elementary news iw ith the value ¢ given in Eq. 25 for A(or the v; value 
given in Eq. 37 for K), the maximum moment in a beam acted upon by a con- 


7 

_ This value may be compared with the exact value obtained from Eq. 7and Table : 
- 3 ‘The variation of the error induced by using the approximate value of Eq. 

43, isgiven in Fig. 5. 


Ina similar manner, using the elementary theory the following values 3 are 
found for maximum deflection and pressure under the loa 


PA _ 58) vy 563 


0. 1.813 


» 


Ne 


12. “Evaluation of Cositicients of Subgrade Reaction,” by Karl Terzaghi, Géotec 
e, Vol, V, 1955, pp. 297-326. 

ee “Uber Die Berechnung der Fundament-deformation,” by Fr. Vogt, Avhandlinger, 
No, 2 Norske _Videnskaps - ~Akademi, | Me ath .- Klasse, 


i 
For 
= 

Eq. 40 yields = ra 

> 

— 
4 

| 
and 
— 


LASTIC SOLID 
_ These values may be compared to the exact value obtained from Eqs. 5 and 9 
with numerical values from Table 2. The errors induced by using approximate 
values of Eqs. 44 and 45 are shown in Fig. 5. 
om Examining Fig. 5 and bearing in mind that the “exact” values of integrals in in 


‘Tables 1 and 2 were computed with an accuracy of + 2%, the following conclu-_ 


4. ‘re beams of infinite length on an elastic- -isotropic semi-infinite sub- 
_ grade are analyzed by means of the elementary theory based on Winkler’s hy- 
B  sgecry with the value given by Eq. , 25 for A (or the wine given in Eq. 37 for 


- ‘FIG, 5. —ERRORS IN THE MAXIMUM BENDING MOMENT M, DEFLECTION Ww, 


dl amount of er error ‘depends on on the e relative ‘stiffness of the beam. ‘Table = 
4 gives the average | and the maximum error P for maximum bending moments, 


OF APPROPRIATE METHOD OF ANALYSIS 


‘The preceding has revealed the fact that the conventional analyete 


is reasonably accurate for beams of infinite length provided a proper value of 
‘coefficient k (Eq. 38) is chosen, However, earlier studies of beams of finite 
length have shown clearly that, depending on the flexibility of the beam, the - 
entional analysis: may also furnish which from 


q 4 
deflections of the beam are underestimated, 
— 
— 


to use the } dimensionless length c characteristics AL in establishing a criterion 
for distinction of cases where the — of the conventional analysis is | 
beams of medium length (0.80 L 6. 


3. short beams (A L < 0.80) | 
at Since the bending moments of a long beam can be computed with an error 


of less than 2% assuming that the length of the beam is infinite, it is not dif- 7 
>. ficult to realize that the statements made on the validity of Winkler’s hypothesis: 
for beams } of infinite length ; are valid without restrictions for long beams too. | 
Also, an analysis of different data available in the literature2,3,4,5,6,7, . 7 
7 clearly shows that Winkler’s hypothesis is not satisfied for short beams. 7 If 

an accurate analysis of such beams is to be made, procedures similar tothose — 

used by De Beer, Ohde, or others must be us‘ used. “(of course a very rough es- 

timate of bending moments in such beams can be made assuming — i- 

gidity of the beam and linear distribution of pressure. eae ae 

pa Concerning beams of medium length (a case of great practical interest din 

to the fact that most of the ange belong to this category) analyses made in- 


TABLE 3. —ERROR IN ANALYSIS OF INFINITE BEAMS ON ELASTIC- motnore: 


a7 SUBGRADES BY THE CONVENTIONAL k- METHOD, 


Influence 


Average error or 


dicate ‘that, although Winkler’ s hypothesis is not strictly satisfied in this range 
of AL, reasonable | accuracy may be expected using the 


if the length « characteristics \ L stays within a certain limit. ayer 
To determine this limit, computations were made of bending moments in 
beams of different ‘stiffnesses Eb I resting on an elastic- ‘isotropic ‘subgrade, 
acted upon by a a concentrated load in the middle. Conventional analysis 
with K-values determined by means of Eq. 38 as well as a more rigorous an- 
eS the subgrade as elastic-isotropic solid were used. Fig. 6 shows — 
the amount of error in maximum bending ‘moment involved by the use of Wink- 
ler’ 's hypothesis. is seen that with increase of characteristics A L,as defined 
_ by Eqs. 36 6 and 38, the error decreases, _ becoming practically negligible for for 
_ Thus, in analyzing a beam of finite length under a concentrated load t by the 
conventional analysis based on Winkler’ 8 hypothesis, reasonable e accuracy is 
‘to be expected if the beam is sufficiently long (AL > 2.25). _ It ca can be shown - 
that the same criterion may be used when several loads are present on the 
beam with the total free length outside of the loads as L. The fundamental fal- _ | 
lacy of the conventional computations of beams on elastic. “subgrade lies in ap-— 


plication of the to short beams. The “superposition meth- 


— 
| 

a 
| 
| 
7 

| 


” for b beams of ties are iil applicable if the “ end- conditioning for forces” 
e ‘small enough compared to the loads acting on the beam. 
These findings have been substantiated by large-scale tests per formed re-_ 
cently at the Georgia Institute of Technology, At! Atlanta, with steel beams resting 


subgrade of compacted silt.. 14 4,150 


_ 1. A rigorous solution of the considered problem indicates that the ratio of 
: pressures and deflections along the beam is practically constant if the beam 
has | an infinite length. The approximate » value of this constant ratio, defined in 


6. —ERROR IN MAXIML MOMENT UNDER A CONCENTR: ATED. LOAD IF THE 
COMPUTATION IS MADE BY THE CONVENTIONAL ANALYSIS | 


the elementary theory of beams on elastic subgrade as the coefficient ati 


- 


grade reaction k, can be computed by Eq. 38. 
When beams: of infinite length resting on an elastic-isotropic half-space 


are analyzed by means. of the elementary theory based on coefficient of sub- 


14 “Beams on Elastic ‘Subgra ade and the WwW inkler’s s Hy ypothesis, A, Vesic, I 
4" ceedings, Fifth Internatl, Conf. on Soil Mech, and Foundation Engrg., Paris, 1961 (i 


1s “Model Studies of Beams Resting. ona Silt Suber ade, B. Vesi¢ and w. H. 
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grade reaction k, bending -moments in the beam a are overestimated, contact 
"pressures and deflections underestimated. The amount of error, depending on 
wn in Fig. 5 and Table 3, 3. —_ | 


> 


«Recommended Procedure 
Size of Beam | Criterion eave. 
| for distinction | Rough Estimates | For Refined Analysis 


Long Beams AL>5.00 Conventional a analysis assuming infinite beam 
and estimating K by means of Eq. 38 


Lé 
Moderately Conventional analysis as as finite beam; K- -value alue 
Beams | 2.25<AL<5.00] by meansofEq.38 = 


Short Beams 0. 80<AL<2, 25 | conventional analysis similar methods 


Short Beams ti treat as De Beer, Ohde, or 
» 4 { 


pigid similar methods 


Concerning beams of finite length, it is shown that. the ‘conventional 
alysis based on the elementary theory is justified if the beam is sufficiently 
long. Depending on the value of the length characteristics AL, the 


given in Table 4 are Fecommended, — 


_ APPENDIX I I.—COMPUTATION OF THE INTEGRAL Jo(0). 


ee to the infinite value of the ir integrand for « “- 0, the integral was divided» 


log 4 


= is, as s before, the abscissa of the point from where on the integrand can be 
approximated bya simpler analytical function representing its upper limit (Fig. 
3).° Both € and M are chosen so as to assure an “accuracy of + 2% in the final 

ve compute Dl” it should be noted that, for ‘small enough values of a, 


1 ae 
c) becomes negligible with av(p). 


2 

TABLE 4,.—RECOMMENDED PROCEDURES FOR ANALYSIS 

4 

| 
Here ¢€ is an abscissa small enough that /() can be approximated by its as- BS ti, 

— 

. 


oy) 50) 


| is integrated using Simpson’ s rule and choosing the number — 


of f points interpolation s so as to assure an accuracy of 2%. 


~ M is chosen great enough ‘so that Pry canbe. disregarded i in com- 


1costM 1 1 é2costM 1 


Evidently for = 0 this expression becom 


Che 


APPENDIX II. —NOTATIONS 


‘The following symbols adopted for use in the paper and for the guidance of 


discussers, conform ny with “American Standard Letter Symbols for 


= width of the beam = = 2b, in inches ~ 


= half-width of the beam, in eal 

- a function of BC c* 1.00 if the distribution of pressure across the 


is uniform; 1. 13 if the deflection across the beam 


Introducing for its asymptotic value, Eq. 47, yields 
J |0.923 - log(2 @)| da .... (49) 
which becomes, by integration, _ 

O] 
4 “parisi 
(51) 
ae 
ag (53) 
q 
x 


fundamental length of the beam = |C Ep |}! 


Young’ ‘modulus for the beam, bunds er square inch» 


= Young’ s for the subgrade, in per square inct 


moment of inertia of the beam, in inches 
J5 


defined in Eqs. 5 through 9 and Eqs. 13 


modulus of subgrade reaction for an infinite beam, in pounds per 


modulus 0 of subgrade reaction fora square plate, in 1 pounds per square : 
coefficient of subgrade reaction, in pounds p per cubic inch 


value of coefficient of reaction for abeam of infinite length, 
a in pounds per cubic inch > 


= value of coefficient of subgrade reaction for a square plate, in pounds 
we length of the beam, in inches 7 
“7 Bending moment in the be beam, in inches 


=a particular value of a ‘_t 


ies a couple acting or on the beam, in pounds rs inch 


a concentrated honk, in pounds 


os pressure at the contact betw een the the beam and the ‘Subgrade, in pounds 7 
shearing force in the beam, in pounds 


= deflection of the beam, in inches 


abscissa along the beam, in inches . 


P= 
increment of along the beams, in inches 


nowy abscissa of the first zero of Jo, in inches | 


Zz any influence in the beam > 


influence in the beam due to a concentrated load P = 1 
sav very sn small value of a 


- slope of the elastic line o of the beam 


/3 

q 
Ko 

— 
— 

“a 

i 

> 

= 

af 


= damping factor- -expression (36). a in in inches~ 
a Parameter sad (23) in 


= a function tabulated f for B > 0. . and given a an asymptotic ex expres-— 
sion for 0.1. 


4 


= Poisson's ratio for the subgrade 
4 
— ae 


aA 
: 
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CONTINUOUS BEAM- -COLUMNS ON ELASTIC FOUNDATION 
Bys. L. Lee,? , ASCE, M. Wang, and J. S. ‘Kao? 


4 


- ‘The analysis of continuous beam- columns on elastic foundation by means 


the slope-deflection equations and related problems are presented. _ The appli 
_ Cation of the moment-distribution method to problems involving large numbers ~ 


of spans is demonstrated by numerical example, and the computation of the 
critical axial force is / illustrated | by the treatment of a particular problem. 


—The letter symbols adopted for use in this paper are defined 

where they first appear, in the illustrations, | or in the text, and a are arranged 


alphabetically, for convenience of reference, in the Appendix, — 
General.— The differential | equation of the curve of a beam- 


< 
on an elastic foundation with uniform cross section? shown in Fig. lis 
= e - 


cox Note,—Discussion open until September 1, 1961. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE, This paper is part : 


of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the © 


American Society of Civil Engineers, Vol. 87, No. EM 2, April, 1961, 
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Evanston, Tl. 


Northwestern Univ., 
Evanston, Ill, 


2 Prof, of Civ. Engrg., 
2 Graduate Student, Civ. Engrg. Dept., Northwestern Univ., 
7 , Graduate Student, Civ, Engrg. Dept., Northwestern Univ., 
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— 
hes. 


in which E 1 is the flexual rigidity, N denotes the axial force, k represents the 7 
- foundation modulus and p(x) is the transverse load, The general solution of 
3 be written inthe form 
= (Ay Ag 2 ™)cos(ax) +(Ag Ag e sin(ax) +y_, s(x). 
_ ich y (x) isa particular integral of Eq. 1 1 and 


a Substituting Eq. 2 in| the 


and 


E 
_leads to four equations, the solution of which yields the constants of integra- 


“tion Aj, A2, Ag, and Aq. ‘Substituting these constants in Eq. 2 and taking the 
ticul 
to the equation of the curve 


N ( 


— 
7 
and 

— 

| 

sin(hax) sin(ax) + y (x) 
=. 
1 
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+l y, 


Sin cos sy afte 
cos (h sin 
sin (hv) ) cos (hv) 
die 


8 °2 


tg = Lo 


and 


lis the of the beam-column, 
Differentiating Eq. 6 with respect to x to x yields the sl slope 


(x) B x) ‘cos (a (ax) + B x) sin(a +Cq 


| Mab! (ea) q 
4 
ft 
x) 
a) 
+ 


Ca 


aC 3+ 


(x) = =C, cos(h Bx) cos(a ‘cos(h ax) sin(e x) 


I 


+ sin(h 8x x) x) ‘sin(h, sin(a 


~ 


a 


= 3 a” + (a> - 3a (16d) 

7 ‘It is obvious from Eqs. 8 that, under any loading condition, once the end 
- moments Mab and Mba and the end deflections Yq and Yb. are determined, the 
equation of the elastic curve is given by Eq. 5. The determination of tl the enc end 


moments and end deflections will be examined subsequently. | 
> 
"Substituting Eq. 6 in the boundary conditions: 


and solving the resulting M Mpa yield the deflection 


equations 
5 <a Studies 1 (Slope Deflection Method),” by G. A, Maney, Univ. of Minn- 
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— 
‘tas. 
and 
13) 

7 1 = (g = Cy ? 2 a B C4 ee ( 14a) 
— 
and Cy = (g2 = a“) C4 2 a B Cy ee © @ ( 14d) 
(x) = Cy cos(h x) cos(a x) + Co 

sin(h 8x) cos(ax)+Cqg sin(hBx) Sinta (xp... 
in which VG pforceand 
p) C4 (a 1 ee ee ( 16b) ‘he 

| — oe 


3 K 


| 


2 


ae 2, 
me + [2 6 sin( 


(0) Ky, (1) -K1 (0)... (20a 


fixed- end moments, and 2 are if acting, on 


- ¢ cos $sin(h 


sin 


for analyzing continuous fi frames, 
(Fig. 1 1) is hinged, that = 0, (a) becomes 


6 “Analysis of Continuous Frames by Distributing Fixed- End Moments,” by H. Cross, 


Mab = K 0, + Sp K + S3 K>~ - Sy K + Mpay -- - - (18a) 

= So K 0. + S; K 0..+S,K b > Ms. . 

In Eq. 18, Mab and Mba are positiyv © ill aClilig ClLOCAWISE OF) UIC member, 6 
and are positive if clockwise and downward, Yq and Yb are positive (Fig. 1). 
— 
= 85 K 85K Mpyy CM ppg (22) a 


@ sin(h 2 v) +z sin(2 


ni The procedure for computing the end deflections and end moments by the 


use of the slope- -deflection equations is best oan by means of numerical 
example, 


shown in Fig. 2(a). 
‘ag support, the > spring constant nt being 30 E vis. 


= 2. 8728 K 6.5862 K - 14, 478 x 107 


Mpa = 9% 3010 K 6, - 21.064 K — 14,.478x 10 K 


inw which K= E Setting x x = 1 in Eq. yields" 


5 
_ For span be, d= 3. 0266 and v= 2,8566 wher eas y_(x) = 0. In this case se Ba. 


Mic = = 5. 5.6684 K K +17 x ( 
setting x = 0 in Eq. leads to 
= 


‘Substituting Eqs. 25, , 26, Pa ‘and 28 in Eqs. 29 leads to 


969 3.7920 14, 478x107 = 0... 
3.7920 + 295.75 - 119,51 x 107 
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inwhich 

| 

ConSiaer tne continuous beam-columns on elastic foundation loaded as 
a 
— 
in 1S Ca 
...(25>) 

8) 


0.86765 x 10 and yp = 


values in 25 and 27 gives 


= Mie = = 1,8688x 10° K ....... ++ 


Finally, substituting the end end deflection of each span in n Eq. 6 


ab (*) = 10" 1 [- 0s Ay x) cos(ay x) /1,0132 cos(h By? x) 


0919 sin(h By cos(a9 x) + + 0.90525 sin(h 8 x) sinla x) +1. wall 


(x) =10- 1 [0.3924 39249 x) cos (ag x x x) sin(ag x) 
= 
0, 39583 3 sin(h Bg +x) cos( (02 x)- 0. 08534 -sin(h Ba, x) sin x x 
in which a, = 2. 3186/1, By = 2. 1510/1, a 0266 / and Bo = 2. 8566/1. 
"The expressions for y "(x), M(x) and V(x), if may be obtained 


- either by means of Eqs, 11, 13, and 15 or by successive differentiations sof 


‘MOMENT - _DISTRIBUTION METHOD 
_ When the number of spans increases, the number of unknown displacement 
- components and equilibrium equations, such as Eqs. 30, increases accordingly. 
To avoid the solution of large number of simultaneous equations in such cases, 


it can be seen from Eq. a. 18(a) that the stiffness of a is K 
. that is, Sy K is the value of Map required tocause 4, = 1 and Op = Ya = Yb = 0, 


at beam- column shown in Fig. 2(a) treated in Example 1 will 1 now be com-_ 
puted by the moment-distribution method, This is done in two stages shown 
in Figs. 2(b) and (c), respectively. . The bending stiffnesses of members ba 
and be are 9.3010 k and 5,6684 k, respectively, as computed in oo. 25 and 27, 


Hence the distribution factors are re 


3010 K + 5.6684 K ye 
6684 K 37866 ae 
“be 3010 K + 5.6684 K 
The carry- ieiaanae — b to a is 0. 30885 as given by Eq. 21 for = 2. 3186 


and = 2.1510. Since end c is hinged, the -over b toc c is is of 


94x 107" 1. 
q 
4 
— 
| 
= 
Ne ena moments Can pe OV OL = 
fined by bq. 41, For a beam-column with the far end hinged, the bending 
of the near endis Ss Kas indicatedinEq.22 
o e near end is 55 Kasi a q. 
nerica 
a 


In Fig. 2(b), an auxiliary force Ry is introduced at ‘point bso a as to re 
to 
strain point b from deflecting. The corresponding moment distribution oper 
ation is shown in Table 1, The distribution and carry-over factors are entered 


TOW The fixed- “end moments due to the uniformly 


in row value of the force at bis 


RL = (0) - = ( 14.769 +70, 1632) (10° 85.401x10°° .. (33) 


in which Vie (0) and Vay p(1) are computed by Eq. 15, 


allow for the of the spring support, point is 
flected downward a distance as shown in Fig. 2(c). The 
sponding fixed-end moments are 


=. 21, 064 x 10°? k 


and are entered in row 5 of Table ... MFa b® and MFba are computed from Eq. — 
by setting = = “Ya =0 and | Mpa b= = MFba = 0, and MFpa is 
mined ina similar manner by Eq. _ 2. The distribution and wey" -over i. 


‘sponding value o of is 


BEAM 

| 
— 
— 


BEAM COLUMNS 
_ The end moments and Vp of the system ‘shown in Fig. 2(a) can now be 
termined by the the ‘Superposition ofthe two —* shown in Figs. 2(b) and (c), 


TABLE i. TABLE FOR EXAMPLE 2 
| 


217,256 

21.064 


pie: Let M, M’ and M” ’ denote, ee the enc end moments of lie: 
three systems, It is evident that they a are related by 


in which | aisa factor determined by the equilibrium condition 


which yields a =0, 39314, Hence, the deflection of point b 


a )- 0,39314 x 1 

and the enc end mo moments are 


17.256 + +a(- 5. 9,559 x 10 107° (40a) 
= K) = 1,8727x 10" K. .(40b) 
_ These values are i n good ag agreement with the the results of Example 1 : 


“CRITICAL AXIAL FORCE 
beam- on elastic foundation.—Consider again the continuous beam- 
columns shown in Fig, 2(a ). In the present instance, however, we 0 and the 
loading consists of Nj No = =0.5 pEI/I? and N3 = I in 
"which pis a parameter, ‘The two equilibrium equations by ‘Eqs. , 29 
apply in this case with the exception, of course, that Mpa: ~ ab(1), and 
(0) are changed. _ Expressing these end moments. and e end shearing forces 


in terms of #, and Yb by means of iss 18(b), Eq. 22 and Eq, 15, respectively, : 


ru 
= R (37) 
als 107 1) = 85.401x10°" —— - a (187.23 . (38) 
— 
4 
5 


a 


substituting 1 the e resulting ex in Eqs. 29 lead to 


(0) 


in 


cos (h vy) sin oy 


+ 


1 


_ 
a 


' 
? © , 


In Bas, 42, the subscript 1 ont: the superscript prime denote items for span 
b, 2 and double primes refer to items =o be ‘be and, by Eqs. 3 and oie 


gene uations ( (Eqs. 

is zero and the — 


— 
= 365 9 81) + - 3 61 10 — 
| cos (h vy) cos + (vp - 3 vy) £2 8 
etka ( - 39, So +t Si) sin(h sin 
= y 2) + (8 3 b tio | ....42) 
x(p) = 2) S4 id S3 £14 + (93 - 3 - 
| 

vy = 5 - Te 
__ For arbitrary values ofthe parameter 
41) are satisfied only when 6, = yp = 0, 


‘ 


_ BEAM COLUMNS 
‘continuous beam-columns remain straight. ‘However, Eqs. 41 have non- zero 
Solutions if the determinant of the vanishes, ‘that | is, 


only unknown in this The of p 
Ea. 44 defines the critical axial force desired. This value can be determined 
7 approximately by the intersection of the D-p curve defined by Eq. 44 with the © 


D= 0, the curve e being plotted by at least three in the 
Thus, 


The intersection of a smooth curve passing through these ’ points ‘with 
D=0 yields the critical value Po = 
=. 


Cantilever _beam- on elastic | foundation such as the ¢ one shown 


free 


FIG. 3. —CANTIL EVER BEAM-COLUMN ON ELASTIC FOUNDATION | 
= 


‘integration AL “Ag, As and Ag in 2 determined by a 
ene 
inthe boundary conditions y(0) = Yar y = = = 0 and = 0. 
the thes se Cc 


n(h ) sin(a x) + (x). 
in which Vn (x) is a by Eq y Eqs. is. 8(a) and on respec-_ 


18 yp). (46a) 


| 
4 
» 
a 
a 
a (82,5) = 344.20 
— 
— 
4 
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_ (y? - ) sin 2 sin? 


_2 2 92 sin(h2 (y2 2) sin? 


- sin(2 9) 
£15 


7 
°15 


re 


‘The expressions for y "(x), y” (x), and y (x) are similar to those given — 

y Eqs. 11, 13, and 15, and may be obtained by replacing the subscripts 2 and 
by 5and6, respectively, 

The slope-deflection equation for a cantilever beam-column on elastic: 


oundation is obtained by setting y’(0) = in Eq. 45 yields 


22 
aC 2 ‘sin(2 0) 


66 | 
q 
— 
....47) 
— 
= (y -39 ) cos(h cos +20 (0 
_ 29 (9 - 37 ) cos (h cos 9 - 2a (y 
 sin(h 2 - Sin an | 
sin 
io 
>| 
| 


BEAM COLUMNS 


2) si sin(h sin sin - 2 ») co cos (49¢) 
[v co cos (h nv) ‘sin 9 + +¢ ‘sin(h wv) ‘cos (498) 


an 


tre ated 


DISCONTINUOUS LOADS 


The saatienine solution defined by Eq.7 can be conveniently used in dealing 
with continuous loads extending over the entire span, In the treatment of dis- 
continuous and y (x) may be ‘Tepresented by trigono- 
)= n 


“Substituting | Eq. _ in Eq. 1 leads to the particular integral 


1+k1° 


7 For a prneenemee load P located at a distance s from the left _ for wm | 


stance, 


“Analysis of Continuous Beams by Fourier Series,” by S. L. Lee, "Proceedings, a 
ASCE, Vol. 83, No. EM 4, October, 1957; Transactions, ASCE, Vol. 124, me at 164. 


Afansactions, 

= 


4 
ae 

Continuous beam-columns whose end spans are cantilevered Can 

the manner these ex) 
| 
a 
se Although Eq. 51 converges rapidly, the speed of converge th 4 
| 

— = 
‘ 


52 leads to the more com- 
pact form of the particular solution 


In involving discontinuous loads, therefore, it is necessary to 
replace Eq. 7 by Eq. 55. 


K n= 


CONCLUSIONS: 

In the preceding, 
which N is a tensile axial force, the sign of the second term on the 
left side of Eq. 1 becomes negative, and the general solution of the resulting 
differential equation is obtained by interchanging o and £ in Eq. 2, a and 8 ? 
being still defined by Eqs. 3, Therefore, all the expressions designated by — 
numerals derived in the foregoing may be used in this case by simply inter- 

_ changing ¢@ and v with the exceptions of Eqs. 51, 52, and 55. In the last three 
expressions, the” signs of the terms containing N’ or p should be changed. 
ma Another analogous problem is that of a beam with uniformly distributed — 
elastic restraint against both deflection and rotation. If the rotational re- . 
straint against both deflection and rotation. If the rotational restraint is pro. 
sos toy (x), that is, the restraining moment in an element dx is 


| 

‘the differential equation” is obta: 


longitudinal pa of unit width of the cylinder, N the axial compressive fence 
acting on the strip, E I should be replaced by ———— and k by E h/r In 7 


the latter, p s ratio while and denote the thickness and the 


“whereas ring ‘stiffeners are the spring ‘support of Examples 1 and 
2. Cylinders with stepped wall thicknesses can be dealt with by | sateen 
the individual sections of constant thicknesses as separate spans, a 
With: reference to ) Examples" 1 and 2, the conclusion may be 
the continuous beam-columns are composed of two or three spans, a solution : 
a may be conveniently obtained by applying the slope- deflection equations di- 


rectly. Whenthe number of s spans increases, however, the moment- distribution 


—_ 8 “Beam or on Elastic Foundation,” ” by M. Hétenyni, The oe, of ‘Michigan Press, Ann 
Mich., 1946, p. 139, | 


4 i 
= 
&g 
a 
_ _ The present analysis Can also be extended to the problem of the axisym- — 
™ j. metrical bending of cylindrical shells subjected to radial and axial loads. In © ane 
| 
&§ 


method is s preferred since the successive approximation technique the 


number of simultaneous e equations needed to achieve 4 solution. ie 


APPENDIX,—NOTATION 


‘The are adopted foruse in inthe paper and for the 


= constants ¢ of integration; 
= subscripts denoting points; 
= Fourier” coefficients defined by Ea. 53; 

= undetermined coefficients; 
Co, C Ce, Cg = constants | of integration by Eqs. and 46; 


= carry-over factor defined by Eq. 21; 


flexural rigidity; 


= thickness of cylindrical shells; ir 


= = = foundation modulus; : 


= length of beam-column; 


= ‘moments; 


end moments at a and b in member ab; : - 


Mpay = fixed-¢ end moments defined Eqs. 20 and 49(d); 


icentrated load; 


oad function; 


reactions at point | al 


= distribution factors at joint bs 


= Cer 


distance of ‘Pf from left end; 


end shearing forces at a and b in ab; 7 


= 


= shear function; 


distributed load; 


= 


— 
or, ; 
4 
a 
— 
4 
k 
&§ 
p(x) 
“a d by Eqs. 19, 23, 49(a) and Fh 
a 
— : 
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= = end deflections of member ab; 


= deflection function; = 


— = deflection of simply supported member ‘when N=0 
and k = - 0 defined by Eq. 52 Eq. 553 


a = particular integral defined by Eq. 7 or Eq. , 55; 


= defined by Eqs. 
constant; 


constants in deflection functions defined by Eas. 


coefficients which are functions of 


= = end rotations; 


= Poisson’s ratio: 


= constants defined by Eq. 49(c), (e), and (f); and 

= 

land £1, 


=a 


— 
— 
(x) 
a8 
2,0, 
os 

= 
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‘ANALYSIS OF FRAMES WITH NON-LINEAR BEHAVIOR® 


H.-S. ANG, M. ASCE. —It is reassuring to know that 
Wilson’ s independent results check ‘ closely the method presented by the writer. 7 
Wilson points out the effect of vertical loads on the primary forces as a 
result of large lateral displacements. This point was purposely omitted in the 
7 “original paper because theoretically, the moment-curvature relationship is a 
- function of the ratio of the magnitude of vertical load to the critical buckling 
load of a given member. 11 The presence of vertical loads, therefore, is only _ 
-— considered when the effect of geometry alone is taken into account, 
However, if the vertical | loads are small compared to the critical buckling 
loads of the columns, the moment- curvature eens for pure bending can 
be used without appreciable error. 
was suggested in the paper that if simple between 
“ments and the end slopes (for 


can be rome, then the total | a for each member can be expressed as in 
equations of continuity « can be formed for each joint: Solution of 
S approach seems to be along 

the same general line. He assumes a otic moment-curvature relation- 

_ Ship, as given in Eq. 12, and obtains moment end-slope relationships given by 

moment-curvature relationships with marked “discontinuity, Eq. 12 
does not give an accurate representation immediately following initial yielding. 

_~Eq. 12 with a=3 and n= =7 gives” too high a resisting moment for a given cur-_ 


--vature in this region as “may be seen in Fig. 13. The moment-curvature rela- 
_ tionships presented as Figs. 3and4arefor pure bending. These agree closely | 
_ with experimentally determined relationships, of members under bending and 


ee: for regions in which measured moment-curvature measurements are. 
available. 12,13 


ao. 


10 asst. Prof, of Civ. Engrg., Univ. of Ilinois, Urbana, Ill. 
a 11 «Static and Dynamic Load-Deflection Tests of Steel Structures,” by F.L. _ Howland, _ > 
Ww, Egger, R. J. Mayerjak, and R. J. Munz, Univ, of Illinois Civ. Engrg. Studies, +, Struc- 
tural Research Series No, 92, February, 1955, 
4Commentary on Plastic Design in Steel: General Provision and Experimental 
— ation,” Progress Report No, 2 of the Joint WRC-ASCE Committee on Plasticity © 
Related to Design, Proceedings, ASCE, Vol. 85, No. 3, July, 1959. eRe boa 
_ 13 «Shear Deflection of Wide-Flange Steel Beams in the Plastic Range,” by W.J. Hall _ 
and N. M. Newmark, Transactions, ASCE, 122,1957, 0 
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oe With the. method proposed by the writer, it is not weaitiihi to determine. 
directly or in a single computation the deformation corresponding to a cm 
load. However, in certain problems related to structural dynamics, it is the 
resistance of a structure’ corresponding to given deformations that is of in- 
terest. It is in this type of problem that the method has greatest application. 
One advantage of computing for the load corresponding to a given deformation 
is that in most ‘Stress-strain: or load- ‘deformation 1 relationships f for “rt 


yield point, “Thus, the strain 0 or curvature to a given stress or 
+. moment is not uniquely defined in n such regions. Such difficulties are not en- 
countered with the writer’s approach. 


ose = 
=i 


— 
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WIND STRESSES 1 IN 


Discussion by Gerrit H. Tosbes 


mation concerning “stresses due to anti- symmetrical wind pressures ina 
hemispherical dome elastically built intoa cylinder.” The result would ina 
TF opengl “their use by the designer.” The solutionis found “on the conservative | 
side” and because the accuracy of the solution “improves very rapidly with 
decreasing thickness ratio, the solution is acceptable inall cases of practical — 


in 


i The Static Loading of Shells. —It appears that the previously cited 


_anti- symmetrical distribution of the load. Although S. -Timoshenko3 
adopted such a sinusoidally varying wind loading, a look at the potential flow 
ry around a circular cylinder, together | with the application of the Bernoulli © 
_ equation would show that the wind load at the sides of the cylinder should be 
far from zero, In fact, the suction at the arene parallel to wind direction 


aims are not quite fulfilled. The presented method of analysis is based on an : 


“edge « of the cylinder, An inspection of the streamline | picture for the flow ‘past a 
a sphere > would afford a similar conclusion. The top of the dome pictured in 

‘Fig. would thus be subjected to a large suction loading. 

‘The least that should be done to aid designers of large shells is to . select 

a realistic loading based on wind- tunnel or water-tunnel tests. The paper by | 
Placido Cicalal? has recognized the uncertainty in regard to the loading. : 
Aside from adopting amore realistic loading from French regulations Cicola’ a 
_ development of an analysis involving atruly general type of loading is a perti- 

nent and an inherently fruitful approach to the problem, 
In 1958 a series of six papers by a Committee on Wind Forces appeared in 
the Journal of the Structural Division. Based on Swiss experiments, loading 

distributions were presented!® relevant to the discussed geometry. This 

_ formation has been compared the authors’ _anti-symmetrical loading in 
‘Figs. 8 through 11. The discrepancy is evident 

_ The Dynamic Wind Loading of Structures.—Even though building codes» 

7 point the way in translating wind forces to static loads the designer should 
remain aware of the dynamic aspects of the problem. Utilizing static loads» 


even if f they correspond to winds of cyclonic or hurricane strength, may give 


4 4 October 1960, by P. Gondikas and M. G., Salvadori (Proc. Paper 2616). lo = 


16 Asst. Prof. of Hydr, Engrg., Purdue Univ., Lafayette, Ind. | a 
ad “Membrane Stresses in Hyperboloid Shells of Revolution,” > by P, Cicala, Proceed- 
an ‘—" “Forces on Enclosed Structures,” by T. W. Singell, Proceedings, ASCE, Vol. 84, 
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—PRESSURE DISTRIBUTION 
wy OVER SPHERE (REFER- 


17, FIGURE 2- 24) 
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PRESSURE. DISTRIBUTION 


PRESSURE DISTRIBUTION 
OVER CYLINDER AND 
HEMI SPHERE UTILISED 


IN F REFERENCE a, pal 


POTENT! 


DISCUSSION 
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‘Atlantic: coast near New York. Structures of novel of relative. 
low mass and of larger dimension than previously contemplated should not be | 
designed without some consideration of the turbulence characteristics of of the 
wind or of the wake flow created by the structure. __ 
> — Certainly this question constitutes a field of interest in itself. _ Despite con- — 
siderable research efforts in many quarters19,20 the picture is not cache 
for even the well studied | case of a circular cylinder. This is not too 
: surprising. A mathematically reasonable description of turbulence must be a 


_ Statistical one, so that ‘securing experimental data, descriptive of | a turbulent 


oy 
11. —EXPECTED DY NAMIC PRESSURE 


TRIBUTION OVER CYLINDER IN TUR=— 
BULENCE FREE WIND FOR 


there exists a lack of basic data on “structured turbulence. 
theless, it may be useful to pursue the question when the adoption of a cor- © 
‘rectly distributed static (mean) windload might be unjustified. For the purpose 
of discussion four situations are distinguished below. In actuality these may a 


not be found separately. For the structures considered, it is assumed that | 
-:19 “The Hydroelastic Behaviour of Flat Plates as Influenced by Trailing Edge Geo ag 
metry,” by A. T. Ippen, G. H. and P,S, . Eagleson, Hydrodynamics Lab., Mz T., 
20 “Fluctuating Lift and Drag Acting on a Cylinder in a Flow at Supercritical Reynolds 
Numbers,” by Pang, Journal of the Aerospace Sciences, Vol. 27, N 


4 a false sense of security. Perhaps the latest example that may be cited is the 4 a 
4 
4 
— | 
7 


‘light material, compositions of large pz panels, tall steel structures, “and so on, 
Structures Placed in a Turbulent Shear Flow.—At distances of more than 

100 ft to 300 ft from the earth’s surface the wind will resemble a homogene- 

- ously turbulent flow. Below this height at which turbulent mixing is enon 
or the ground the flow will be a turbulent shear flow. The mixing process 
-_(gustiness”) will be intermittant (“squalliness”), but still entirely random. 
Both the scale relative to the structure andthe intensity of the turbulence oa 
a 3 be known for a rational selection of design loads. Unfortunately, too little 
= been published of the work done inthe field of micro-meteorology to per- 
‘such selection. Insight into the problem may be gained from a paper21 

“= Robert H. Sherlock, F. ASCE. Analytical techniques to treat random loading - 

are available. To mention a typical example: the case of structural panels ~ 
located in a flow with random disturbances has been treated by J. Houbolt. 22° 
‘In view of the ‘random 1 character of the loading itis probably ‘adequate to 

atop an equivalent static load. It should be realized, however, that the static 
7 suggested by building codes” may need to be modified to account for 
differences eographically) in turbulence intensity or for the ratio of 
_ scale of maliadanae to that of the proposed structure. The recommenda- 

- tion? by R. H. Sherlock, -F. ASCE, isa ‘step il in the correct direction. It is not 

felt that the adoption of a judiciously chosen static design loading would be 

imherently unsatisfactory, 
Structures: Placed in a Wake Flow. _—Not all is random > 

‘in character. The wake behind blunt obstacles, whether this is a skyscraper 

or a ridge in the form of a long factory building will contain a certain amount — 

of periodicity as well as turbulence of alarger scale than present in “natural” a 

_ wind. Structures placed in such wake flows are said to be buffeted. The po- 
“font danger is of course that a predominant frequency component in the 
"pressure spectrum of a wake flc flow would coincide with a frequency of 
Information concerning the pressure spectra | of wake flows in addition to 
e~ knowledge of the natural frequencies of the ‘proposed structure may be © 
"necessary for the design» of stacks, shells, towers, and SO on, located close 
to other structures. Analytically the problem would not be too involved23 once ; 
sufficient information is available concerning the wake flow in question. 
Fortunately it has been found24 that going down stream the periodic com- 
_ ponent in wake flows. degenerates into random turbulence quite e rapidly. Except — 
in cases of closely | placed objects, such as a row of smoke stacks, it is felt 7 
that the danger of resonant buffeting will remain small, Further research 
are needed, however, to substantiate the statement. 

_ Structures Creating Marked Wake Flows.—Of greater “consequence may 

be the paneen lift force, FL, exerted on structures due to their own wake | 


22 “= “On the Response of Panels Subject to a Flow Field | Containing Random Distur-_ 

bances, ” by J.C. Shock and Vibration Bulletin Part I, Office of the U.S. Secy. 


AWA 


_ 2 “Wind Forces on Structures on Structures: Neture of the Wind,” by R.H. H, Sherlock, — 
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Gi a 24 “On the Development of Turbulent Wakes from Vortex Streets,” by A. Roshko, — 


"DISCUSSION. 
flows. _ The “usual, nee. for serious wind induced vibration is that von 
Karman vortices are shed periodically at a frequency given by the Strouhal 


ys 


in which Sis the Strouhal number, a object; Ws 
refers to the Strouhal frequency of vortex shedding, in radians per second 
D denotes the characteristic transverse dimension of the structure, in feet; 
and U is the free stream velocity, in feet per second. . When the frequency Wg 

coincides with a natural frequency, of the structure, resonance would 

ome Neither the magnitude of S nor the magnitude of the lift force, Fu is 

From a fluid mechanics point of view afew remarks can be made that may 
“bring the problem into better focus. Th. von Karman’s, Hon. M. ASCE,25° 
analysis of the stability of a two-dimensional vortex array neither dealt with 
the mechanism of vortex shedding nor furnishes” an adequate description of 
wakes for Reynolds numbers larger than (in the case of circular cylinders) 
R= U D/v = 300 in which vp is the kinematic viscosity of the fluid. It is a 

_ matter of record, however, that the mechanism of vortex shedding persists 

up to the velocities at which the boundary we gd becomes seprinennge prior to 

- separation (at the “critical reynolds number”). Therefore, even though the 
wake becomes increasingly suppressed by general tur- 


of the wake for subcritical Ps. en adopts the relations as a convenient and ~ 
__ description of the wake structure it is readily shown’ that: ‘lel 
s cp = constant 


D= -—», - = drag | coefficient = 
D : 
"Experimental one for a circular. cylinder verify reasonably the relation : 
given in Eq. 67 with a value for the « constant of 0.16 to 0.22. For non- circular : 
- eylinders26 the constant appears to vary from 0.15 to 0.30. Eq. 67 may serve. 
q 7 asa simple guide in considering geometric modification of a two- dimensional 
blunt cylindrical structure in order to avoid periodic lift forces of an unde- -— 


sirable frequency. Many drag data are either available or at small 


As mentioned previously, with increasing Reynolds number the wake periodicity : 
7 is increasingly suppressed by general turbulence leading to an increasingly 
random dynamic loading. The increase of general turbulence in the wake should 
associate with a break-up | of the two-dimensional wake structure that exists ; 
7 at low flow velocities. The scale of the turbulence relative to the structure 


= is de sing R. What remains of the discreet vortex fila- 
’ ments along the trailing edge will be generated in a mutually random phase 7 


25 “Uber den Mechanismus des Flusigkeitsund Luftwiderstandes,” by Th. von Karman 
and, H, Rubach, Phys. Zeit. No, 13, 1912, pp. 49-59. 


“Low Speed Drag of Cylinders of Various Shape,” Delany and N, E,Soren-— 
‘NACA Technical Not Note 3038, 
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“a relationship. Consequently, the net ‘lift force will decrease with increasing 


"aa ‘It thus appears that in judging the effect of a dynamic wind loading associ- 
ated with the wake flow behind a structure, both the geometry of the structure © 
_ and the scale of the wake turbulence constitute independent variables that 


- must be considered. Unfortunately, little is known about the wake flow struc- 


7 ture as a function of the Reynolds number to assess adequately the role of 


a static wind load would appear acceptable as long as a non-periodic, tur- 
_-bulent wake may | be anticipated. This must be qualified, however. A A mechanical 
_ linear system may well be excited to vibration at its natural frequency by a 

random force. Only a small part of this force does effective work, But, if the. 
- force is large, this may yet lead to large stresses. , Because lift forces in- 
Be oni with U2, winds of cyclonic strength n may y conceivably impart ‘sufficient — 

kinetic energy to a structure to destroy it. This possibility depends on the 
_ elastic and internal damping characteristics of a structure. 
Aeroelastic Vibrations of Structures Creating Marked Wake Flows. —Inthe | 
- discussion so far the tacit assumption was made that the elastic response does 
4 not influence | the flow field. However, it is only by considering thi this response 
that a consistent description can be for the cases of really 

"serious wind- or water-induced vibrations. 

Detailed ‘measurements and analyses ot ‘resonant response of flat 
plates in water tunnel18 | permitted the conclusion that, dependent on the 

‘geometry between the separating boundary layers, the effect of large amplitude 
a motion (that necessarily occur at natural frequencies of the structure) is to ; 

restore a degree of two- ~ dimensionality in the wake. The | elastic response thus 
"provides the aero-or hydroelastic effect of phasing the broken up vortex fila-_ 
a ments. The net lift force increases, giving rise to larger amplitudes of vibra- | 
tion and again an increase in lift, The higher the velocity, U, the 
_ more dangerous this self-excitation becomes because F o U2, The a ell 

‘motion will govern the frequency of vortex formation over a limited to a wide © 
range of free stream velocities mainly depending on the structural geometry © 


near the trailingedge 
_ These observations are made having in mind subcritical R’s and essentially . 


= o-dimensional structures” such | as the author’s cylinder. 
_ dimensional structures, such ; as the hemi-spherical dome proper treated by i 
the authors 7 dependency of the wake structure on the structural response ~ 
q = be assumed to exist. Factual knowledge about the wake dynamics of — 
= dimensional bodies limited. is known, however, that the base 
4 


“pressures: just behind the trailing nies are higher than for two-dimensional 
7 configurations, From A. Roshko’s interference tests23 it may be deduced that 
1 ¥ low base pressures are essential to the formation of strong vortices. Hence, 
_ a is tentatively concluded that for three-dimensional structures the chances 

on aero- or hydroelastic vibrations is less than for t two-dimensional geome- 


tries. For the number region (that is different from 


tm 
_ lation of lift forces along the trailing edge does not preclude the intermittant _ ie 
occurance of large lift forces. But their duration usually is too short to cause 
tm 
; 
i 
: 


larger structural motions ar are likely to effect the turbulence, too. It probably 
arrests the trend to a smaller and smaller scale for the wake turbulence. If 
this happens, afurther increase in wind velocity will bring | increased excitation 

7 because Fy, o U2 and despite the fact that F, remains a random force. The 
= situation thus has the prerequisites for destructive Self- excitation ‘wer - 
“consistent with known facts. It ‘makes 2 a ‘difference ‘that structural materials 
are used, Larger deflections and lower internal damping make steel structures 

_ more prone to serious vibration than concrete structures. The natural fre- _ 

“ quencies of a a structure can be as important a consideration as the computed 
static stress distribution. The permissible magnitude of deflections under 

_ loading may depend on considerations dealing with turbulent flow. = 

It is realized that the preceding remarks were not directed towards the | 

7 analytical core of the authors’ paper. Rather, the claim as to the immediate | 


practical value of the analysis has been 
L 


1. Analyses directed towards” the determination of wind stresses ina | 
_ structure ‘should at least be based on experimentally obtained distributions — 
of the mean aerodynamic pressure force. Fluid mechanics can not furnish 
entirely rational descriptions of the pressure distribution over two- or three- _ 
iv 2. Analyses directed towards the determination of wind stresses | ina 
_ Structure should attempt to incorporate as general a loading as is practical. a 
i Even for one e and the same structural geometry the loading pattern may vary 
“surface irregularities or proximity of other structures. 
‘engineer should always realize that the translation of turbulent 
7 wind pressures to a static design load may constitute a dangerous simplifica- 
- tion. . That the real problem is a dynamic one may be of little consequence in 
some cases, but canbe determinative inothers, 
4. Questions that require attention in the case ofimportant or novel struc- 


7 tures for which wind loading is animportant consideration, include the follow- 


a. The turbulence characteristics of the oncoming air r stream; — 


velocity, generated in the wake of the structure; 
The ratio, as of wind velocity, of the scale of structure 
structural the internal damping 1g and the 
_ magnitude of the elastic r response during | cers motions at natural fre- 


_ These questions ‘pelong to the more difficult o ones in the field of structural and 


fluid ‘mechanics . Future structural ‘designs will no doubt require answers, 
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BEARIN 
ae BEARING 1G CAPACITY OF FLOATING ICE SHEETS® 


on by 1 L. W. Gold, R. H. Wood, N. C. Lind, D. T. Wright 


w. GOLD.5- 5. —In 1957 a that arose the discussions that took 

; place at the symposium sponsored by the National Research Council of er i 

a on the bearing capacity ¢ of ice was” a ‘survey among the pulp and paper com- * 

- panies of Canada on their field experience with the use of ice covers for stor-_ 3 
: logs. It is common practice in this industry to build up the thickness of a 
_ the ice cover of selected lakes. and rivers, and when this thickness is adequate, 

: to transport the cut tim iber onto the cover by various types of trucks 3 and 
sleighs. Many thousands of cords of wood are placed onto the ice in this way. 
‘Individual loads can be as large : as 70 tons to 90 tons. _ Such extensive us use of | 
ice covers is certainly a rich source of field experience. 

- _ When Meyerhof prepared his paper, the records were available of ice 
failures obtained from the survey during the first winter. Some of these Ob- | 
‘servations are shown in Fig. 8) along with observations on maximum re 
normal loads. Since that time, further information has been obtained, part of 
_which has been published. 6 The additional observations on ice failures fall = so 

into the same region as those given by the author, — eo 

Co It was recognized early that many of the failures occurred for ice e thick- 

nesses that would normally have been considered adequate. Greater attention 

. 7 was then paid to the possible cause of the ice failure for each breakthrough. 


To date most of the breakthroughs that could not be attributed to thin ice were © 
‘considered due to either (1) excessive vehicle speed, (2) repetitive use of the | 
_ ~Same road on the cover, or (3) the presence in the cover of a water layer be- — 
With: respect to the first point, it is the ex) experience of men men responsible for 
_ traffic on ice landings that an empty truck travelling at a speed between 15 
- mph and 20 mph will “often cause more audible cracking of the cover than a 
7 fully loaded truck travelling at speeds not exceeding 5 mph. With respect to : 
the second point, it is the practice of some companies to change the location 
of a regularly used road approximately every 2 hr. It is considered that this” 
action has been beneficial in reducing the “fatigue” effect of the traffic. 
Relative to the congee! of loaded vehicles that pulp and paper companies 
place covers, a breakthrough is not a common occurrence. During 
the winter of 1960- 61, a daily record was kept at four large ice landings of 
the number and size of the loads _— onto the ice, along with the ice thick- | 


October 1960, by G. G. -Meyerhof (Proc. ‘Paper 2627). 


5 Head, Snow an and me Sect., Div. of f Building | Research, Natl, Research ’h Council, Ottawa, 7 
PA. “Field Study on the Load Bearing Capacity of Ice Cov ers,” by L. W. Gold, Woodl: ands 
;, Pul and Paper Magazine of Canada, Vol. 61 No. 5, 1960, 11-153 — 
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ness. Only a ‘part of are available as of March 1961. Those 


that” were received showed that the average loads would lie between the lines" 
fy = 50 and fp = 100 in Fig. 8. The maximum loads would extend into the region 
_ between fh = 100 and fy = 200. The weights of the loads, including the truck, | 
were between 10 tons and 75 tons and the traffic density at times exceeded 150 
trucks in 24 hr. No ice failure’ has ian senna as yet in this survey of the 
dt In the analysis. of the cheorvations diated from the survey on ice failures, 
» $ no attention has been paid to the influence of geometry of loading on the bear- . 
ing capacity. It was considered that it would be unrealistic to take into ac- 
_ count the load geometry in view of the lack of knowledge o on the : ice e properties — 
"and the errors associated with estimating the load and measuring the ice 
a thickness, The rectangular space defined by the wheels of a truck has an area - 
of between 100 sq ft and 200 sq ft. ‘This area is equivalent to a circle with 
_ radius» “a” equal to about 7.5 ft. If the ice has a thickness between 25 in. and 
30 in., the ratio of this radius to the characteristic length L would be about 
(0. 15. For the majority of vehicles and ice t thicknesses normally e encountered, — 
equivalent a/L computed in this way would lie in the range 0,12 a/L 
The | author « does state: that the usual wheel spacing f for a unit has no in 
fluence | on the e collapse load in practice. It would be appreciated if the author 
would enlarge on this opinion and in particular comment on whether it is suf- 
ficient under n normal conditions to estimate the collapse load from an equiva-_ 
“Aen value for -a/L obtained by letting ‘ “a’ ’ equal the radius of a circle enclosing 
the same area as defined by the wheels of the vehicle. 
- Some time ago, the writer made observations on the loads required to 


cause radial and circumferential cracking in thin ice sheets formed in a tank | 
The loads were applied at the center of the sheet over an area 0.5 in. in ~ 
; ad radius. The edges of the plates were frozento the side of the tank. A hole was — 
_ drilled in the sheet near the edge sothat the water under the ice was not con- 
7 d fined. The loads were applied fairly rapidly. For most of the tests the cracks © 
; formed within five minutes after loading was begun. The weight of water that 
: came up through the drilled hole and onto the top of the e plate was ; always less 
“than about 2% of the applied load. 
a ‘Table 1 are. given the thickness of the plates, their characteristic 
— length L, the load required to form the radial and circumferential cracks and 
the radius of the circumferential crack. It was assumed that the ice was in ~ 
the elastic state for each test. Eqs. 1 and 38 were used to compute the tensile | 
stress, A for the radial and the tensile stress fh for the. circumferential | 
_ erack situation respectively. The ratios of the radius of the contact area of 
the load to the characteristic length of ice were in tecadl range 0.013 < a/L << 
; 0. 033, and the ratios of the radius of the tank to the characteristic length were 
inthe range 
It is realized that the values obtained for S, and fh are approximate only, : 
; _and likely on the high side due to the finite size of the plate. The characteristic dj ; 
L was computed from the expression L = 40n3/4 obtained from field 


: sonia hole being less than half the radius of the circumferential crack. In 
test number 6, the load was held at 17 Ib for 15 min and then increased until il 
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: the circumferential cracks formed at 36 lb. For tests 7, 8, and 9 the load re- 7 
quired to produce radial cracking was maintained until failure. For test num- > 
. 7 this was 100 min, test number 8, 40 min and for test number 9 neither > 
the circumferential cracking nor oe had occurred when the oad w was s 


Thickness, in Py, in} S,, in 

in inches inches pounds | pounds per pounds 
square inch| | " per square 


_in inches 

square inch} inches |wedge angle, 


pounds per sector of 


o 


= For two t tests, the plate was cut inhalf and one of the halves into quarters, _ 
7 The load was : applied at the edge in. the center of the half plate and at the tip 
a’ of the 90° sector plates. The results are given in Table 2. Eq. 38 was used to q 
- compute fp. For the last test in Table 2, the first radial crack occurred when 
the load was 49 lb. Using Eq. 14 gives for Sy a value of 250 psi. -dicaepriadigll 
The crack pattern for the edge loaded plate had the shape described by 
Meyerhof. The ratio of the distance along the edge from the load to the crack 


__piscussioN 
test 
0.56 | 25.7 | 375] 242 | 955 | 293 | 10 
5 | 0,52 | | 283] 207 
7 28.2 | 30.9) 159 | 
291 
0.73 31.6 m/2 180 
0,73 (31.6 | 4 
Average 8 
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to the same ‘distance along the bisector | was between 1.4 and 1.8. This com-— 
pares” quite well to the same ratio for the contours of constant deflection — 


ne tables it is seen that the caine for Sy are about —- 
s 


ame as those obtained for fp. All the ee are within the range snonpendl 


ports the validity ‘of the “expressions by the author for estimating 
% the collapse load of floating ice sheets under short time loading conditions 
--The writer was on a lake on the first day that loads of wood were being 
placed onto. the ice. Cracking could be heard near the truck as the trucks 
travelled over the ice at less than 5 mph. This cracking was probably radial 
only, no evidence of circumferential cracking was observed. The total load 
was about 40 tons and the ice thickness about 25 in. For the condition that no 
_ cracking occurs and assuming L to be 450 in. and a/L = 0.20, the value com- 
puted for S; using Eq. 1 is about 365 psi and for fp using Eq. 38 is about 120 
. psi. If the / cracking that did occur was radial, the ultimate strength of the ice 
in the cover would, — therefore, be between | these two values and this is in 
agreement with the laboratory rer 


H. WOOD, 8 A.M.1.C.E.—The author’s proposals to investigate the load 
carrying capacity of large ice sheets by means of rigid-plastic theories of g 
plate” action is interesting and no doubt intended to be a first approximation. a 
In point of fact the actual collapse loads are likely to be somewhere in between 
the predictions of ey loads according to plastic and elastic theories 
respectively, perhaps» an interaction diagram could be plotted, after the 
style used by W. Merchant, M. ASCE,? 10 for predicting the 
collapse loads of building frames, if only 1 more e test results were available. 

_ As discussed elsewhere by the writer!1 one of the main problems with 
plastic design of. slabs” is to discover a satisfactory yield criterion on 
he base the mathematical treatises. Itis perhaps surprising that Mr. aol 


Johansen ‘criterion is much a much more extensive 
of known solutions available. 12 In addition, the work of E. H. Mans- 
‘fieldl3 | is appropriate to | the use o of the Johansen criterion, in 1 predicti 
worst patterns of collapse by variational methods, while the writer!! has 
also found lower- bound solutions for collapse modes similar to that nine ; 


_, Dept. of Scientific and Industrial 
9 | “Frame Instability in the Plastic Range,” by W. Merchant, Symposium on Plastic >i 
Theory of Structures; Cambridge, British Welding Journal, Vol. 3, 1956, p. 366, _ 
“The Stability of Tall Buildings,” by R. H. Wood, Proceedings, E., Lond 
Vol, 11, September, 1958, p pp. 69-102; also Vol, 12 . April, 1959, pp. 502-522, yy a 
Il «plastic and Elastic Design of Slabs and Plates,” by R.H. Wood, Thames and Hud y. 
“Brudlinieteorier,” by K. Ww, Johansen, -Teknish Forlag, Kopenhagen, 1952, 
oun “Studies in Collapse Analysis of Rigid- o-< Plates with a Square Yield Diagr “am,” 
by E E, H, Mansfield, Soc., A., Vol. 241, 1957, pp. 311-338, 
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(2 +F)Mo based 

on Tresca’s s criterion. In the extreme case, thereis the well known difference 


that a clamped plate should collapse under a point load when P = 2 7 Mo by 


Tresca’ s criterion, but when P = 4 Mo by Johansen’s criterion. Such differ-— 


It would appear doubtful whether the yield moment could be based on fp h2 /4, 


‘that a rwx equal yield strength and equal ductility in tension and compres- 
A. J. Ockleston14 and the writer!) have found considerable evidence of 


Johansen load, by t cietee a no- in in- concrete theory into a kind of self- 
prestressed ‘slab in which cracking is inhibited. There is then an instability © 
type of failure, but the highest moment that can be generated 
. : ‘by membrane compression is of the order of fy h2/12 . This suggests that one 
a ‘should be careful before adopting a criterion of maximum-plasticity-: in- eal 
theory h2 /4) in order to interpret the test results. 


Both membrane action and elastic foundations have one thing in common— -, 


‘ments, That being the case ‘es author’ s inclusion of the upward water pres- 


_ sure is interesting but unnecessary. It was a foregone conclusion with rigid- 


_ plastic theories that the same starting collapse load would be given as for the 
” ‘standard proofs using the Tresca criterion, because atthe start the deflection 
- must be zero by plastic theories. As regards the evidence from elastic theory 


the writer would appreciate a helpful definition of the “characteristic” length al 


This interesting treatise by the meee, toe together with information a 
that is- accumulating from other sources, seems to point to the — 
that there is. considerable scope in the future for theories of failure 


‘that do not ‘depend on a well defined yield point. 


15 . 16_ 
UN. LIND, and D. D. ‘T. WRIGHT. —It worthwhile to under 


the states of loading considered in the paper. Ii is recnacrt that the —— 
moments» in the ice sheet be developed at reasonably small deflections, and in 
- any case before the state of loading is noticeably modified through pmol 
: gence of the ice sheet. Although this requirement is most important for corner 


‘in which. ceili of the ice from other causes is likely ‘to occur, The follow- 
ing example serves to establish some approximate limits for the cases in 
_ which the deflections should be taken into account. 
—_ _ Consider a strip load of width 2a = = 0.4L on the edge of a semi- infinite ice | 


sheet. - Assuming elastic action, Eqs. 18 and 19 give the maximum deflection 


14 «| “Load Tests on a Three Story Reinforced Concrete Building in Johannesbur, g,” 7 
by A. J. Ockleston, The Structural Engineer, Vo No, October, ; also Vo 
k J. Ockl Th 1 Engi , Vol, 23, No. 10, b 1955; also Vol, 
24, No, 10, October, 1956. Cann 
"15 Assoc, Prof, of Civ. Engrg., Univ. of Waterloo, Waterloo, Ontario, Canada, | 


16 een of Civ, Engrg., Univ. of Waterloo, Waterloo, Ontario, Canada, a 
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a gence in the active zone of the ice sheet is estimated at one fifth to one half : 
of the thickness h of the ice sheet. Using Wo < h/3 | as a criterion for negligible 
deflections, Eq. 1 yields h > 10 in. for relatively weak ice (fp = (100 psi ) to d 
h > 21 ft for strong ice (fp = 500 psi). ‘Consequently, neglecting the influence 
of deflections in the elastic analysis for short- -time lo loadings must be viewed ed 
with caution, especially for ice of high strength. 
en _ For long-duration loads, using a typical value of fh = = - 30 pi psi, the preceding» 
criterion indicates that the elastic deflections | may be neglected in practice 
(h > 1 in). However, the inelastic deformations would presumably be many, 7 
4 times the deflections indicated by Eq. 1. In the absence of data regarding 
a reasonable first estimate would be approxi- 
mately three times the elastic deformations . Then, with fh = (3) (30) = 100 psi, 
_ Eq. 1 yields h > 10 in. This result indicates that deformations associated | 
the | inelastic behaviour of ice are worthy of further study. 
An unresolved matter of importance is the strength of natural ice. Analyses : 
of behavior will induce little confidence without reliable knowledge of the — 
ranges but not always consistently. influence of factors such as 
a. density, and solarity could be determined quite readily through load : 
tests on small samples in the f field.17 Unfortunately, there is a marked scale 
effect in the strength | of ice, | and small sample tests cannot now be readily 
yg related to the strength of natural ice sheets. A theory for scale effect on ice 
_ Strength thus appears to be sorely needed. V. V. Lavrov18 has proposed ‘such 
theory, with hypotheses that are not entirely convincing. 
_ Rational determination of the load capacity of floating ice sheets is s thus 
seen to depend on gorse for determining strength | in the field, a ee 


7 “Recommended Standards for Small-Scale Ice Tests; ” by T. 


"actions Engrg. Inst. of Canada, Vol, 2, No.3,p.112,1958. 
: 5 “The Nature of The Scale Effect in Ice and The Strength of The Ice Sheet,” by = 
‘vz. Lavrov, Soviet Physics Doklady, Amer, Inst. of Physics Vol. 3, No. 5, 934, 1958, 


=. strength of the material. As the author indicated, strength varies over a 
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6 and L2 =v E' he / v1 2k yields q | 
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ELAST ASTIC BYD PROC -EDURES® 


Discussion by Daniel Frederick — 


DANIEL FREDERICK, 6A AM. ASCE.- —The writer wishes to © point out that 
the method used by the author for the elasto-plastic analysis of structures is _ 
practical and time- -saving -when compared with several other methods that 
a might have been used. ‘Dean T. Mook? considered the elasto-plastic analysis 
of bents by the complementary energy method taking a rectangular cross 
section and a material having a bilinear stress strain diagram with equal | 
values in tension and compression and all other conditions the same as those 
“used by the author. The complementary energy unit t length along the beam 


for the elasto-plastic cross section, 

was » integrates along the bent from one elastic zone to pany This with ~—_ 

usual expression for complementary energy in the elastic portions gave the 

total , complementary energy. The usual methods based on complementary : 

energy were then used to derive a set of non-linear algebraic equations in- 


: volving the maximum value of a for each plastic zone. Although the method 
> is exact, the ‘solution | of these equations i is tedious | and not readily amenable 


equations are not fully developed. Other methods, such as using transformed © 
sections lead to the same difficulties. Based on this experience, it appears 
that the approximate method used by | the author is one of the best for the 


sve 


analysis of bents also. 


= December 1960, by Annabel Tong (Proc, Paper 2690). 
prof, of Engrg. Mech. Virginia Polytechnic Inst., g, Va 
ae T« Analysis of a Bent Constructed from a Strain- ~Hardening Material, ” by Dean : 
_ Mook, Thesis submitted in partial fulfillment of the requirements for the degree of — 
Master of Science at the Virginia Polytechnic Inst, at Blacksburg, Va., in1960, 
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